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ABSTRACT 

1. The constant A, derived from Cowling’s formula for apsidal motion in a close 
binary, is computed for 19 systems. Many of the values are rough, but all are small. 
The largest value found from reliable data is o.0o18—very much nearer to a Roche model 
(0.000) than to a homogeneous one (0.750). Refinements of calculation would diminish 
the calculated values. The values of K certainly differ from star to star. 

2. The “photometric ellipticity” derived from the variation outside eclipse depends 
on the darkening at the limb and the variation of surface brightness with gravity, as 
well as on the geometrical form. The gravity effect approximately doubles the observed 
variation (eq. [14]). The gravity effect can be well studied by this means; the geo- 
metrical shape and internal constitution cannot be found in this way. 

3. The light-curve resulting from an eclipse of two prolate, uniformly bright, el- 
lipsoidal stars of different shapes can be very closely represented by the eclipse of two 
stars of the same shape but of different dimensions. The outstanding residuals simulate 
the effect of a small degree of darkening at the limb. It is impracticable to find the 
difference in shape from the light-curve. 

4. It is recommended that the ellipticity of figure calculated on dynamical prin- 
ciples from the Roche model (or perhaps with K = 0.01) should be used in determining 
the geometrical elements of eclipsing variables. 


I. APSIDAL MOTIONS AND CENTRAL CONDENSATION 
1. The theory of apsidal motion in close binaries appears to have 
been cleared up by the recent work of Cowling’ and Sterne.? Walter 
now agrees with their conclusion that for gaseous, and therefore 
almost perfectly fluid, bodies, the figure at any moment will be very 
nearly that of equilibrium under the instantaneous tidal forces. 


*M.N., 98, 734, 1938. 


?M.N., 99, 451, 1939. 
3 Zs. f. Ap., 17, 320, 1939. For solid bodies, Walter’s former conclusion that the apse 
regresses is correct. 
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The periastron advances. The principal term in its mean motion, 


n’, is 

n’ m, R,\$ mM, 

= (1416) +x, (2) (1 + 6%), (1) 
where 1 is the orbital mean motion, R, and R, the radii of the stars, 
r that of the orbit, m, and m, the masses, and K, and K, constants, 
depending on the density distribution, which have the value 0.750 
for homogeneous stars and 0.000 for the ‘‘Roche model” with all the 
mass at the center. Terms of higher order in R/r and in the eccen- 
tricity increase the motion. If the mass ratio, radii, and apsidal mo- 
tion are known and we set K, = K, in (1), we obtain a mean value of 
K which gives information regarding the internal constitution. 

2. Of the observational data, m,/m, is rarely directly determina- 
ble. Unless R, and R, are decidedly different, a rough value will suf- 
fice, and this may be obtained from the relative brightness (assum- 
ing, for example, that m varies as L‘’). 

Errors in R, and R, are much more serious. Owing to observa- 
tional selection (see below), the systems in which apsidal motion has 
been established show primary and secondary eclipses of comparable 
depth, rarely exceeding o™6. These are just the conditions under 
which the determination of the relative radii of the stars from the 
photometric data alone is practically indeterminate, even with a 
good light-curve. The ratio of the surface brightness is determined, 
and also the sum of the radii, from the duration of the eclipses (so 
long as these are not too shallow). 

As a typical case, we may take the depths of minima equal, and 
the sum of the radii as 2R, and set R, = R(1 + a) and R, = R(1 — a) 
We then have 
(1 + a) 2/3 


I-a 


and we find, for the coefficient C of K.R’/r5 in equation (1): 


a 0.0 0.2 0.3 0.4 
1.00 1.50 1.86 2.43 
34.0 a5 05 40.0 56.7 
(C/34)'/S. . 1.000 1.009 1.032 1.068 1.108 
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The last line gives the factor by which R must be multiplied to give 
the same value of K for equal stars. A solution assuming equal radii 
will give the maximum value of K. This will be a good approxima- 
tion unless the radii are decidedly unequal. 

3. Unless the light-curve is really poor, the duration of the eclipse 
at half the maximum depth can be found from it with an uncertainty 
of a few per cent. The whole duration, between the unobservable 
outer contacts, must be calculated from this. If 6, and 4, are the 
longitudes (assuming a circular orbit) corresponding to the first con- 
tact and half the maximum loss of light, we have‘ 


sin? 0, 
sin? 4; 


= x(k, Ao ’ 


where & is the ratio of the radii and a, measures the maximum 
eclipse. This function has different values for uniform disks, for 
darkened disks with the small star behind, and for the same with the 
small star in front.5 Call them x, xo, and x’. If 1 — d, is the depth 
of the minimum with the small star behind, and 1 — X, that of the 
other, we have for uniform disks 


Assuming various values of a,, & and x. may be found. For darkened 
disks we have 


ai‘Q(k, 1) = aiQ(k, as) . 
With any assumed value of aj, the values of k, ao’, xo, and x,’ may 
be found with the aid of tables.‘ 

For typical cases in which the minima are of equal depth we find 
the results given in Table 1. 

For uniform disks all the solutions give nearly the same light- 
curve and the same computed duration of eclipse. When the eclipse 
approaches internal tangency of the disks, x. becomes smaller, which 
diminishes the computed length and counteracts the effect discussed 


4 Russell, Ap. J., 35, 326, 1912. 5 Ap. J., 36, 392-394, 1912. 
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above. In this case, however, the eclipse chord may be considerably 
less than the sum of the diameters of the disks, the computed radii 
will be greater, and K will be smaller. 

For darkened disks, if the stars are unequal, the durations of the 
two minima at half-depth will be decidedly different, though the 
orbit is circular and the extreme widths are the same. The average 
value of x, for the two minima is less than for uniform disks for the 
deepest eclipses, but is greater for the shallower ones. 


TABLE 1 


HALF-WIDTHS OF ECLIPSES 


UNIFORM DARKENED 
(1—Ay) 
ao k Xo Ge k xo 
0.50 I.00 | 1.00 6.15 I.00 | 1.00 1.00 5.44 5-44 
0.40 0.80 | 1.00 4.43 0.80 | 1.00 0.800 4.64 4.64 
0.90 | 0.895 4.47 ©.90 | 0.865 | 0.840 | 4.43 4.81 
1.00 | 0.817 4.43 1.00 | 0.789 | 1.000 3.80 5.00 
0.30 0.60 | 1.000 3.76 0.60 | 1.000 | 0.600 4.24 4.24 
0.70 | 0.866 3.74 0.70 | 0.856 | 0.625 4°82 4.34 
0.80 | 0.755 | 3.73 | 0.80 | 0.766 | 0.693 | 3.98] 4.39 
0.90 | 0.707 | 3.69 | 0.90 | 0.687 | 0.795 | 3-74] 4.40 
1.00 | 0.655 3.54 1.00 | 0.601 | 1.000 ce 4.44 
0.20 0.40 | 1.000 3.30 0.40 | 1.000 | 0.400 3.94 3.94 
0.60 | 0.707 2.31 0.60 | 0.720 | 0.486 3.76 4.05 
0.80 | 0.57 3.26 0.80 | 0.586 | 0.656 2" Rie 4.05 
0.90 | 0.535 3.22 0.90 | 0.535 | 0.772 3-34 4.02 
1.00 | 0.500 3.00 I.00 | 0.467 | 1.000 2.82 3.98 


It follows from this survey that in almost all the cases in which 
apsidal motions have thus far been observed the assumption of uni- 
form disks of equal size will lead to a maximum value of the constant 
K in equation (1). 

In systems of a “typical Algol” character, with a large faint com- | 
panion of small mass, almost the whole of the apsidal motion arises 
from the ellipticity of the latter, and more detailed knowledge of the 
elements is necessary. 

4. The observational determination of the apsidal motion is more 
delicate. It is never safe to assume that a periodic term in the time 
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of primary minimum is due to this cause, for light-equation arising 
from motion of the close pair in a large orbit can produce almost ex- 
actly the same effect. Moreover, many eclipsing pairs show more 
complicated irregularities, which are obviously not simply periodic 
and are thus far unexplained.® Only when the secondary minimum is 
observable and shows displacements in the opposite sense to the 
primary can apsidal motion safely be inferred. The absence of stars 
with shallow secondaries from our present lists is thus natural. 

In a few cases a complete revolution of the apsides, or more, has 
been observed, and the motion is well determined. When both the 
interval and the widths of the two minima may be well determined 
from the observations, and also the rate of change of the interval, a 
direct determination of n’ may be made, as in the case of HV 7408.7 
Considerable new information is available for other cases, which are 
therefore discussed in detail.* 

From the interval of minima we have 


us 


= — {ecosw(1 — wo} . 


If only the times of minima are available, we must wait until the 
curvature of the plot of A against ¢ becomes sensible before n’ can be 
found; but it is possible to set a minimum for n’, namely, the value 
which gives the greatest curvature consistent with the observations. 

a) For example, the Harvard observations of RU Monocerotis? 
(P = 345847) give the intervals between minima as shown in Ta- 
ble 2. The residuals O — C are the departures of w from uniform 
motion at the rate given below. For e = 0.40 they are inadmissi- 
bly great; they are tolerable for e = 0.45 or greater. A period of 600 
years (n’ = 0°6) is about the minimum permissible and will give the 
greatest permissible value of K. 

6R. S. Dugan and Frances W. Wright, Princeton Contr., No. 19 (in press). 

7M. B. Shapley, unpublished. 

8 An independent discussion of some of these has been made by Sterne, M.N., 99, 
662, 1939. 


9 Shapley, Proc. Amer. Acad. Arts and Sciences, 66, 469, 1931 (Harvard Reprint 
No. 72). 
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When (as in this case) cos 7 is small compared with R, the dura- 
tions of the two eclipses are approximately in the ratio (1 + e sin w)/ 
(1 — esinw), or 1.7 for e sin w = 0.26. Though, as Shapley says, “‘the 
widths of the minima cannot be estimated with sufficient precision to 
give a value of e sin w,”’ his light-curve is hardly consistent with so 
great a difference as this. For unequal darkened disks with the small 


TABLE 2 
APSIDAL MOTION IN RU MONOCEROTIS 
€=0.40 €=0.45 €=0.50 
|e cosw 

w O-C w O-C w O-C 
BOO ware 2469 | 0.394 | 10°80 | —2°6 29°80 | | 38°0 | —0°%6 
2.65 297 19.6 | +2.2 33.2 | +0.5 41.1 | +0.4 
2.60 27.4 | +3.0] 38.0] +1.5 44.8 | +1.2 


star in front at the eclipse nearer periastron the ratio might be re- 
duced by about 12 per cent (Table 1), but it would still be at the lim- 
it of the admissible, and the period of 600 years may be adopted. 
With so large an eccentricity spectroscopic observations would 
settle the question—and are much to be desired. 
b) Miss Jones’s careful study of V523 Sagittarii’ gives the follow- 
ing normals: 


Date cos w O-C, O-C; 
+0.o11 —0.004 0.000 0.000 
— .0o60 + .004 + .oo1 + .oo1 
— .107 + .oo1 — .003 — ,003 
—0.143 —0.006 —0.001 +0.003 


From these, with various assumptions regarding the eccentricity, 
we find 
(1) = 0.15; w = go° + 2°10 (¢ — 1902.5); period = 171 years 
(2) 20 1.41 (t — 1902.0) 255 
(3) 0.25 1.12 — 1901.0) 322 
'© Harvard Bull., No. gog, 10, 1938. 
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The first period appears to be a little too short; the other two give an 
excellent fit. 

Good mean light-curves are available for both minima. The one 
which comes late is narrower, showing that the periastron is advanc- 
ing. A uniform solution with equal stars, and e sinw = 0.15, givesa 
satisfactory representation of both, and the elements R, = R, = 
0.222, 1 = 85°0, L, = L, = 0.50. This curve should correspond 
closely to the mean value e cos w = —o.08, and hence to e = 0.17. 

Any considerable increase in e sin w gives a greater difference in the 
computed widths of the minima than the observations permit. The 
values e = 0.17, P’ = 200 years, may therefore be adopted with 
some confidence—though for unequal darkened disks e = 0.25 might 
be permissible. 

c) For YY Sagittarii we find, in the same manner, 


Date cos w O-C: O-C, O-C; 
—0.047 —0.004 —0.002 ©.000 
.072 + .004 + .oo1 + .oo1 
.103 + .007 + .002 + .oo1 


(1) € = 0.15;w = 129091 + 1°27 (¢— 1917.8); P’ = 285 years 
(2) .20 116.5 + 0.785 460 
(3) 0.30 107.2 + 0.476 757 


The observers state that ‘“‘no certain difference is now found in the 
widths of the two minima, though about 1880 when w was go’, there 
was probably a measurable difference.’ This is hardly compatible 
with e = 0.20, which makes the ratio of the widths 1.5. We may 
adopt e = 0.17, P’ = 350 years. 

d) Of the two stars observed by Dugan and Miss Wright,” V526 
Sagittarii shows a change in e cos w from —o.18 to +0.14 in 7oooP, 
or 36.7 years. With e = 0.20, 0.25, 0.30, P’ = 125, 162, 205 years. 
The ratio of widths of the minima is about 1.3 according to these ob- 
servations, and 1.4 according to O’Connell."? We may adopt e = 
0.23, P’ = 150 years. 


™ Shapley and Miss Swope, Harvard Bull., No. goo, p. 9. 


'2 Princeton Contr., No. 19 (in press). 13 Riverview Pub., 1, 13, 1935. 
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For V346 Centauri the range in e cos w is —o.12 to +0.09 in 
2250P = 38.9 years. With e = 0.15, 0.20, P’ = 160, 223 years. 
From the rather discordant observations the minima do not differ 
much in width, so that we may adopt the first. 

e) The observations of GN Normae show that e cos w = +0.16 + 
0.0018/ (¢ in years). The widths of the minima are given by Kruyt- 
bosch" as 0.15P and 0.14P, so that e sin w appears to be small and 
negative. The data are rough and partly inconsistent but may be 
represented by e = 0.25, P’ = 680 years. 

f) Finally, there are two systems which show a displacement of 
the secondary, without sensible change. For SS Lacertae (P = 
14°42) the motion may be expected to be very slow (sec. 8). KT 
Centauri (P = 4713) shows e cos w = +0.20—changing by prob- 
ably less than 0.02 in 40 years. According to Uitterdijk,‘’ one mini- 
mum is 0.82 as wide as the other, which would make e sin w 0.1 and 
P’ > 1200 years. But this difference in width is no more than could 
be produced by limb darkening of unequal stars with e sin w = o. 
In this case, w might have ranged from 335° to 25°, while e cos w re- 
mained between 0.19 and 0.21, which would give P’ > 290 years. 

g) For RZ Cassiopeiae a small periodic displacement of the pri- 
mary minimum, of period 18 years, has been detected by De Sitter." 
This corresponds to an eccentricity 0.013, too small to be detected 
spectroscopically; and the secondary is too shallow to determine 
whether it oscillates in opposite phase. This case has been included 
for completeness. 

In a few instances it is possible to determine the apsidal motion 
from radial velocities. Luyten'’ has shown that most of the motions 
previously announced are illusory. Advance is probable for a Vir- 
ginis and 13 Ceti A. In 29 Canis Majoris w was the same in 1916 and 
1935, but the possibility of two complete revolutions in the 19 years 
cannot be excluded. The periastron of AO Cassiopeiae appears to 
make half a revolution (or perhaps 13 or 23) in 8.1 years. 

5. The data, as finally collected for 18 systems, are summarized 


B.A.N., 6, 233, 1932: NOUS AGN 
5 B.A.N., 6, 247, 1932. 17 Pub. Minnesota Obs., 2, 29, 1935. 
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in Table 3. They are listed approximately in the order of reliability, 
with the best at the top. The values of m.,/m, given in parentheses 
have been estimated from the relative brightness; those in square 
brackets are rough guesses—-whose error, however, will not seriously 
affect the results (see sec. 2). The densities p, and p, are in grams per 
cubic centimeter, P’ is the apsidal and P the orbital period, and C is 
the coefficient of A in equation (1). The references in the last col- 
umn give the source of the orbital elements. 

6. For the last four stars the apsidal motion has been determined 
by Luyten from spectroscopic data.'"* The determination of R/r is 
good for 29 Canis Majoris, poor for AO Cassiopeiae, and very poor 
for a Virginis. For 13 Ceti A there are no eclipses, but approximate 
values of these ratios may be obtained indirectly. The close pair is 
a component of a visual binary. Let the masses of the principal com- 
ponent and the spectroscopic and the visual companions be p, ux, 
and Then 


+ 4, + x2) = ps, 


where Y is the visual period in years, a’’ the apparent mean distance, 
p the parallax, D the spectroscopic period in days, r the correspond- 
ing mean distance in solar radii, and C = 1.34 X 10-7. If L is the 
visual luminosity, J the surface brightness, and R the radius, in 
solar units, 

L=JR = Ip", 


where log / = —o.4 (m + 0.15) (m being the star’s visual magnitude, 
and —o.15 the sun’s magnitude for p = 1’’). Eliminating among 
these equations, we find 


R (V\2/3/1 + x, + x,\"3 "1/3 J—1/2 
(5) 


All the factors here are known except J, which may be estimated 
from the effective temperature. 
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For 13 Ceti, a’ = 0.244, Y = 6.91, m, = 5.7, whencel = 0.0046: 
Luyten’s discussion’? gives m, = 0.95, mM, = 0.32,m, = 1.05, whence 
x; = 0.34, x, = 1.1. Then R/r = 0.178J—?. The Harvard spectrum 
is Go, whence J = 1. The spectroscopic companion is probably an 
M dwarf as Luyten suggests, and its radius about half the sun’s. The 
separation of the close pair for mass 1.27 is 7.4 solar radii, and R,/r = 
0.07. These are the tabular values. The Mount Wilson spectrum is 
F7, corresponding to a temperature about 330° higher, on Kuiper’s 
scale, and J = 1.25, which makes R, = 0.167, K = 0.13. 

The apsidal motion is uncertain, owing to unexplained discord- 
ances in the radial velocities. New spectroscopic observations are 
badly needed. 

7. The principal result of this investigation is that the constant K 
is always small. Only in one case is the computed value as much as 
one-tenth as great as for a homogeneous body, and this value is very 
uncertain. Novalue derived from the betterobservations exceeds 0.02. 

This cannot be a result of observational selection, for rapid apsidal 
motions are easier to detect, and these lead to large values of A. The 
approximations made in the discussion (assuming the stars equal, 
and taking the shortest permissible apsidal period) also increase the 
computed K. The more complete theory (which Sterne has applied 
to those cases which warrant it) also leads to a smaller value of K 
than the present first approximation. 

The number of cases now available is clearly great enough to be 
statistically significant. There is, therefore, very strong evidence 
that the density concentration within the stars is high. 

It is also evident that the stars are not all built on the same model. 
Conclusive proof is afforded by GL Carinae and V523 Sagittarii, for 
which the periods and radii are nearly the same, while the apsidal 
periods are 25 and 200 years.” 

8. The value of K corresponding to any given density model may 
be found by integrating Radau’s differential equation by quadra- 
tures. For polytropic spheres Chandrasekhar finds the values given 
in Table 4. The .atio of the central to the mean density from n = 1 


19 Ap. J., 78, 225, 1933- 
20 The same conclusion has been reached in a more detailed discussion by Sterne 
in M.N., 99, 662, 1939 (kindly communicated in manuscript). 
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to n = 4 is given within 5 per cent by the empirical approximation 
p./p = 0.825/K. 

For the stars just mentioned, the values of ” are 2.9 and 3.9, cor- 
responding to 48 and 460 for p./p. Those for the other well-deter- 
mined cases lie between these limits. The lowest value for a poorly 
observed case is 2 = 2.0, p./p = II. 


TABLE 4 
n ° I 2 4 
0.750 0. 260 0.1446 0.0741 0.0144 ©.00134 
I .000 3.290 5.991 II.403 54.18 622.41 
kp-/p 0.750 0.853 0.864 0.845 0.782 0.823 


The two cases of apparently “‘fixed”’ apsides in eclipsing pairs are 
explicable by values of K within the same limits. For SS Lacertae, 
if K = 0.02, the apsidal period comes out 33,000 years. The spec- 
troscopic data are consistent (within their present wide uncertainty) 
with similar values of K. Agreement for AO Cassiopeiae may be se- 
cured by assuming that the periastron has made 13 or 23 revolutions 
instead of 5 revolution. 

Further photometric and spectroscopic observations of these stars 
and of any others which show evidence of eccentric orbits are ur- 
gently needed. Until they have been obtained, the question why 
apsidal motion has been so hard to detect in pairs where the com- 
ponents are nearly in contact must remain unanswered. It may be 
noted, however, that for some of the “‘good”’ systems in Table 3 the 
interval between the nearest points of the stars hardly exceeds their 
diameters. 

II. PHOTOMETRIC AND DYNAMICAL ELLIPTICITY 


It has been maintained that the photometric observations of 
eclipsing pairs outside eclipse show ellipticities which indicate that 
the stars are nearly homogeneous. The relation between this ‘‘photo- 
metric ellipticity” and the actual ellipticity of figure has been dis- 
cussed by Walter” and recently by Luyten,” but something more 
may be added. 

* Schriften der Kénigsberger Gelehrten Gesellschaft, 3, 55-100, 1931. 

22 M.N., 98, 450, 1938. 
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There is general agreement on the forms which are to be expected 
for a pair of stars revolving in a circular orbit and rotating with the 
same period—at least, to the first order in the ellipticity (which 
alone will be considered in this paper). Let m, and m, be the masses; 
take the radius of the orbit as unit, and let a,, 6,, and c, be the semi- 
axes of star 1. Set 


ay, => Cy = ba . (2) 
The ratio of centrifugal force to gravity at the equator is 


m, + m, 3 


Yr mM, I 


Then 


M2 
= a3(1 + 2K;), 


Uy 


(3) 


+m 
= ——— a3(1 + 2K,), 
2 mM, 


where K, is the parameter depending on the internal constitution, 
already defined—and similarly for star 2. 

The light-changes which arise from the rotation of an ellipsoid de- 
pend upon the degree of darkening at the limb and the degree to 
which the surface brightness J varies with the “‘local”’ gravity g.?3 
According to von Zeipel*4 and Chandrasekhar,’5 J should be propor- 
tional to g; but the proof of this depends on the assumption that 
there is no internal circulation of matter; and, in any case, it applies 
to the total radiation. Changes in radiation of a given wave length 
may be greater or less.” We will therefore set, for the apparent sur- 
face brightness of an element of the disk, for which the normal is 
inclined at an angle y to the line of sight 


J = 2+ 4005 y+ 78), (4) 


23 Takeda, Mem. College Sci., Kyoto Imperial U., A, 17, 197, 1934. 
24 M.N., 84, 665, 1924. 

M.N., 93, 573, 1933- 

26 Cf. Kopal, Ap. J., 89, 323, 1939. 
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where g, is the mean value of gravity over the surface. If, as is 
probable, the radiation is nearly of black-body quality, the value of 
y should increase for the shorter wave lengths. The light received 
from the star is, then, 


L= fJ cos ydo (5) 


over the visible half of the ellipsoid. 

If A, B, and C are the amounts of light sent in the directions of the 
axes a, b, and c, then, to the first order in the ellipticities, the light 
sent in any direction will be 


L = C sin? 6 + B cos? 6 cos? t + A cos? 6 sin? ?¢, 


where the direction is toward declination 6 and hour angle ¢, for an 
observer at the end of axis 0. 
When x = y = o (uniform disks), A:B:C = bc:ca:ab; and we 
have 
L, = + u — v+ 0 sin? 6 — u cos? 6 cos? (6) 


When x = 1, y = o (complete darkening),?’ 
L, = + $(u — v) + 8(v sin? 6 — u cos? 6 cos? #)}. (7) 


By Clairaut’s Theorem” the variation in gravity from pole to equa- 
tor on a spheroid is 
Ag 
+ 
+e 
where ¢ is the oblateness. But « = $¢(1 + 2K). It follows that at 
any point on the ellipsoid 


gQ (8) 


27 Ap. J., 36, 397-399, 1912. 

28 Clairaut’s proof of this applied to the case when the distorting force arises from 
rotation and results in a zonal spherical harmonic of the second order; but, as any 
spherical harmonic of this order can be built up from zonal harmonics with different 
axes, the extension of the proof is obvious. 
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where g, is gravity at the distance 7, from the center. To the first 
order, 
y—b 


cos? B sin? — vsin? B, (9) 


where ) and @ are longitude and latitude defined by the equations 
x’ = acosBsind\; y’ = bcos B cos ; = csin 6, 
where x’, y’, 3’ is a point on the surface. 
In the integral for C, ydo is the projection of the element do on 


the plane z’ = 0, which is ab sin B cos B dB dX. 
When x = o and y = 1, we then have 


a J f I u cos? B sin sin? 8) 
sin cos dp 


I 


2K — 
2K +1 


rab], 1+ 4 — 30) 
= 1+ u— —— 


For projection on the plane y’ = o, 


ydo = ac cos? B cos \ dB dd , 


2K ) 
B= wacJ, {1+ (am — 42) | 
sK \ 
= I+ (u = v) 
and, similarly, 
A = rbcJ + ——— (ju — 
We then have 
K 
Ly = (1+; 
10 
+ SETA sint 5 — cos" cost) 
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If we set x = 1, y = 1, in (4), the integrands for A, B, C, must be 
multiplied by the appropriate values of cos y—that is, the direction- 
cosines /, m, n, of the normal. These are proportional to «/a?, y/b’, 
z/c?. To the first order 


1=cosBsindX (1 — u cos? B sin? — v sin? B), 
m = cos cos X (1 + u cos? B sin? — v sin? B), 
n= sin 8B (1 + v + u cos? B sin? A — v sin? B). 


We then find 
2 


L, = 3rb*J, I — v) 
(11) 
\ 
+ pier (v sin? 6 — u cos? 6 cos? t) | 
5 


By (8) the gravitational effect should formally disappear if K = 
and this substitution reduces (10) and (11) to (6) and (7). 
We have, by (4), 


+ y(1 — x) 24 gy -cosy =. 


Hence, in the general case, 
— yL.+ yr — + xyZ,, 


or, substituting from (6)-(10) and neglecting the small constant 
term in u-2, 


2-—K 
— u cos? 6 cos? t) f 
so that, in general, we may write 
r+ ( 2-—K 
+ (LEE) r+ sin 6 


— u cos? 6 cos? t) f ; 


(12) 


Here u and v are the actual ellipticities of the surface. 
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If our ellipsoid is a component of a binary, w and v are (1 + 2K) 
times the values for a Roche ellipsoid, namely, 


of = +m) (13) 
I 7m,’ I I 


We have also 8 = (1/2) — i, t = @ (the orbit longitude), so that 


(1+ + yt 2k 
(14) 


\ 


(v{ cos? i — sin? cos? 6) 


A similar equation holds good for the second component. 
The combined light of the system (outside eclipse) may then be 
written 
L = (L, + L.)(1 — C sin? cos? 6) , 


where we may set L, + L, = 1. Then, to the first order, 


6x; 
C= (: + + y, + 2K, — 99:K,) Liu; 


(15) 
If the light between eclipses is expressed by the usual equation 
L=a+bcos@+ cos? 6, 
we would have 
C sin?7, (16) 


were it not for the effect of “reflection’’ discussed by Eddington”? 
and Milne.*° Harmonic analysis of Milne’s formula for the phase 
variation gives, for the combined light of the components viewed in 
the orbit plane, 


L =1+ (Lia3 — L,a?)(0.347 cos 6 + 0.005 cos 36) 
+ (L,ai + L,a3)(0.073 cos 20 + 0.002 cos 48) , 


29 M.N., 86, 320, 1926. ¢ 39 M.N., 87, 43, 1926. 
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while Eddington’s formula substitutes (0.333 cos 8 + 0.000 cos 34) 
and (0.060 cos 20 + 0.004 cos 48). 
The photometric ellipticity, corrected for reflection according to 
Milne, is then 
C’ = C+0.15(L,a? + L,a?). (17) 


It is preferable to compute this term rather than to derive it from 
the value of 6, as is often done, for two reasons: (1) when LZ, and L, 
are comparable, the reflection term in cos # becomes small, but that 
in cos? 6 does not; (2) there is often an observed asymmetry in the 
light-curve, represented by a term in sin 6, and we cannot be sure 
that the whole of the term 0b cos 6 arises from “reflection.” The 
“Roche ellipticities” u; and wu; may be calculated if the elements and 
masses are known. Equation (15) then gives us one relation with six 
unknown quantities. 

When the components are nearly equal, it is probable that x, y, 
and & are nearly the same for both; but for a typical “Algol system,” 
where the companion is larger, cooler, and less massive than the pri- 
mary, it is probable that the darkening coefficient x differs for the 
two, and possible that the internal-constitution coefficient K differs, 
while nothing is known about the gravitational effect y. 

Even if the parameters are the same for the two stars, we have 
one equation and three unknowns, and it is necessary to find two of 
these from other data. 

From sufficiently good observations of the whole light-curve, the 
degree of darkening, x, may be found, though it is very hard to de- 
rive separate values for the two stars. 

We have seen in Part I that K is small, while there is no doubt 
that x is often large (its theoretical value for total radiation being 
about 0.6), and y may well be greater. Hence, the photometric ellip- 
ticity appears to be quite unsuited for the determination of internal 
constitution. It should, however, in future, give a good determina- 
tion of the gravity effect and so an independent datum for the test- 
ing of astrophysical theories. 

The most recent discussion of the relation between dynamical 
and photometric ellipticity is Luyten’s.** He has taken the gravity 


M.N., 98, 459, 1938. 
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effect into account by adopting a mean value of a correction factor 
but has allowed for darkening by the inadequate method of adopting 
the value of b/a given by each computer, without re-examining his 
assumptions. An exact comparison of his results with the present 
discussion is therefore impracticable. 

Walter? has tabulated the photometric ellipticities on a uniform 
scale and has corrected them for the reflection effect. For uniform 
disks, and without gravity effect, he expresses the ratio of the photo- 
metric to the Roche-model result by a factor K,, which is 3 for this 
and 7.5 for a homogeneous body, so that K, = 3(1 + 2K) on our 
scale. For the mean of 39 systems he finds K, = 6.52,or1 + 2K = 
2.17. On his assumptions this would indicate a low central condensa- 
tion; but the present analysis gives the equation 


6x 
(: + y+ 2K — }yK) = 2.17 


(where the results will be some sort of average values for the two 
stars). On various assumptions this gives 


/ / 


1 3/2 I 1/2 I 
4 ° ° o 1/2 1/2 I 
K=0.58 0.37 0.18 0.38 0.14 —0.08 0.II —0.17 —0.43 0 ° ° 


5 0.36 


y = 


This shows the enormous range of uncertainty. If we set x = 0.6, 
we find 

ner 0.67 — 2K 

For B and A stars the limb darkening should be comparable to that 
for the sun in the infrared, and hence small. As most eclipsing varia- 
bles are of these types, the gravity effect may have nearly the theo- 
retical value. 

The gravity effect, if present alone, would make the center of the 
disk darker than the limb, at the time of eclipse, and so modify the 
light-curve.* From equations (3), (4), (8), and (g) it follows that, 
if we could see star 1 through star 2, the apparent surface brightness, 


32 Schriften der Kiénigsberger Gelehrten Gesellschaft, 3, 65, 1931. 
33 [ owe this suggestion to a conversation with Dr. Kopal. 
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in the absence of the ordinary limb darkening, and neglecting 2, 
would follow the law 


For the reflection effect, at conjunction, Milne’s formula*+ may be 
written in our notation 


The apparent laws of darkening are not the same—nor is either 
one of the usual form; but the ratio of the change from center to 
edge of the disk is 
Lym, 
Lim, 


(3K — 6) a. 

For equal stars and small k, the gravity effect will be somewhat the 
greater of the two, but the resultant of the two will give a limb 
brightening of the order of 0.1 at most, which would be masked by 
the ordinary darkening. The effect would be greater for the large 
faint component of a typical Algol system; but this would affect 
only the curve for the secondary minimum, which is too shallow to 
be determined with great percentage accuracy. It is doubtful, there- 
fore, whether this effect is observable. 


III. DIFFERENCES IN FORM OF THE COMPONENTS 
It has been recognized for forty years*s that, when the compo- 
nents of an eclipsing pair are ellipsoidal but differ in mass and den- 
sity, they will differ also in form. To avoid complications in deter- 
mining the elements, it is usually assumed that the components are 
similar in figure. The effects of this simplification have been dis- 
cussed by Walter,*° who concludes, from consideration of a particular 
case, that the effect on the light-curve is similar to that of darkening 
at the limb but is of opposite sign. 
34 M.N., 87, 50, 1926. 35 G. W. Myers, Ap. J., 7, 8, 1898. 


36 Schriften der Kiénigsberger Gelehrten Gesellschaft, 3, 55-100, 1931 (especially pp. 
82-89). 
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A fuller treatment, however, is desirable and turns out to be 
fairly simple, provided that it is assumed that the two ellipsoids 
present uniform disks and are prolate, with their major axes in the 
line of centers. 

Let a, and a, be the major semi-axes of these ellipsoids, and b, and 
b, their minor semi-axes, a, and b, being the greater. The minor semi- 
axes of the projected ellipses will be 6, and b, and the major d, and d,, 
where 

d? = ai(1 — sin? cos? 6) , 
e, being the eccentricity of the meridian section, 7 the orbital inclina- 
tion, and 6 the longitude from mid-eclipse. The projected distance of 
centers 6 is given by 


6? = cos?i + sin?7 sin? 6, (18) 
whence 
di = bi(1 + ué?) , (19) 
where 
é? 
u = (20) 


and similarly for d,. 
Now let a be the fraction of the area of the smaller disk which is 


obscured by the larger (as for similar disks), and let 


56 = d, + (21) 
If 2s is the common chord and 
sin = sin g, = (22) 
the areas of the segments cut from the ellipses by this chord are 


_A, = bid, — sin COS ¢:) , 


A, b.d,(¢. — sin Y2 cos $2) 


whence 


bid, 
Ta = Ai +. A; = 9, = SiN Cos + — sin cos ¢:) (23) 
b.d, 
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Also, 


p = COS — (1 (24) 


These equations are exact and, with (19) to (22), suffice for the pre- 
cise computation of the light-curve in any given case. 

We are concerned with the question of how closely such a curve 
can be represented by another due to the eclipse of two similar ellip- 
soids, by proper adjustment of the elements of the system. We may 
first find a pair 2 of similar ellipsoids such that the values of 6 at the 
external and internal contacts are the same as for the original pair 1, 
calculate the light-curve for pair 2, and then adjust the elements so 
as to shift this curve into the best practicable agreement with the 
original. We will confine ourselves to terms of the first order in 
ta. 

Then 


d, = b,(1 + 1,6?) . (25) 


Let b{ and b} be the minor semi-axes of pair 2, and w’ their common 
prolateness. Then, by (21), 


5 = + pb. + (ub, + pu.b.)é , (26) 
5’ = bt + p’b, + (u'b + . (27) 
If 6 = 6’ at external contact (p = p’ = 1), we have 


bh —b, + — b, = {b,(u, — u’) + b,(u. — u’)}(b, + 


At internal contact i 


bi — b, — — b.) = {b,(u, — u’) — b,(u, — u’)}(b, — 


whence 


by — by = by + — + — u’)} , (28) 


b, — b, = — + + — w’)}. (29) 
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If, as usual in the theory of eclipsing variables, we set 


k= b.? (30) 
we have 
k’ — k = k(b?2 — + — 2u’). (31) 


If, at any other time, 6 = 6’, (26) and (27) give 
— — — wu’) + (2b, — pb.)(u. — u’)}. (32) 
Let L, be the light of star 1 of pair 1, when seen end-on (6 = o). 
Its light at any other time will be 


b, 


= + u,6?) , 


and that of the system outside eclipse will be 


Lo= L,+ + 


As usual, we may assume L, + L, = 1. 
For pair 2 
L=L+L,+ + 
We will have Li = L, if Li + Li = 1, and 
“= Lu, + Lay. (33) 


The figure to be attributed to pair 2 is thus defined by the “photo- 
metric ellipticity’ derived from the uneclipsed part of the curve. We 
now have 


uy — = — U2) ; u,—u' = —L,(u;— (34) 
During the principal eclipse, of star 2 by star 1, the light will be 
(35) 


To rectify the light-curve so that it is flat between minima, we must 
divide this by L, (or L3), obtaining 


L, =*1 — aL, {1 + (u;, — u’)6?}. 
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The rectified loss of light during totality will be L.{1 + (u, — u’)&} 
and not strictly constant. For pair 2 it will be simply L’. 
If these are to agree at internal contact, we must have 


L, = L,{1 + (u. — u’)(b; — b,)?} ; (36) 
if, on the average, during totality 
L, = + 3(u. — u’)(b, — b,)?}. 
Adopting the first assumption, 
L,=1— abt + (uz, — — — (37) 
For the eclipse of pair 2 we have 
Lr LiA(p’, k’), (38) 


where the function A is defined by the standard equations 


I 
p’ = cos — (1 — cos ¢}) , 


I 
wA = — sin cos + — sin ¢; Cos ¢}) , (39) 


sin = k’ sin 
From (23) and (24) we have 
p = cos — — cos ¢:){1 + — u2)87} , 
Ta = — sin Cos g, + — sin cos ¢;){1 + — u,)87} , 
whence 
Ta = rA(p + (1 — cos — u2)6?, 


I 
+ (er — sin cos (1 — 
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From (39) it follows that 
A 
T = — (er — sin (40) 


whence 


a = tA(p, k) + (yg, — 2 sin + sin ¢; Cos — U2)6?. (41) 


ra’ = rA(p, k) ; sin = p) (es — sin k) (42) 


or, by (30)-(34), 
ma’ = rA(p, k) + — (1 — p?) sin — (2 + pk)L,) 


Combining these results, we find 
L, — Li = — + , (44) 
where 
mC, = —2(2 + pk) sin — + (¢: — sin ¢1) 
(g: — 2 sin + sin Cos (45) 
+ + pk)? — (1 — k)*}a, 
= asin — — 74 *) — sing) 
(y: — 2 sin + Sin ¢; COS . 


The elements of pair 2 which are subject to adjustment may be 
taken as b/ and ’, and the inclination as 7’. 
We have 


62 = cos? t’ cos? 6 + sin’? 6, 


} 
Also, 
Ge 
a 
| 
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where @ is an orbit longitude, which we must regard as an observed 
datum, not subject to change. We then have 


I 


tel 


1 — 6” = cos? 6 sin? i’ ; dj’ = — F 


Hl = 
= + p’)(1 + 


Hence, neglecting w’d6’, etc., 


ds’ = (l' + p’)db, + bidp’ + (48) 
dp’ = —dl (l’ + p’) +p py +P cot z’di’ . 
Hence, by (40), 
da’ = X +sin gi’ + p)¥ + (49) 
where 
X = cot ; Y = ( 6 cot ) 
(50) 
= dl’ J 
or 
cot i’di’ = ; = + b7X) ; 
(51) 
= Ank'(Y + + 


If x, y, and gz are the coefficients of X, Y, and Z in (49), we have, for 
the rectified change in light of pair 2 at any time, 


dL} = —Lida’. 
Equating this to L, — L;, we have 


aX + y¥ + 2Z = —b2(u, — + . (52) 


a 
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These equations are best discussed numerically. Table 5 gives the 
relevant data for two typical ratios of radii: k = 0.8 and k = 0.6. 
By adjustment of X, Y, and Z the computed curve for the new pair 3 
of similar ellipsoids can be made to cut the (hypothetically) observed 
curve for the dissimilar stars at any three desired points. These have 
been chosen to correspond to cos g, = +0.6, 0.0, —o.6. The results 
are as follows, the tabular unit being b7(#, — w.): 


k=0.8 k=0.6 
L2=1 I2=1 
—0.25 +0. 29 —1.12 +1.73 
+ .18 — .21 —o.14 +0.45 
+0.68 —0.83 +1.44 —2.59 


The residuals, in the sense O — C (pair 1 minus pair 3), are given in 
the table. For other values of LZ, and L, the results may be obtained 
by simple interpolation. When LZ, = L., Y, Y, and Z and the resid- 
uals become very small. That is, the light-curve produced by eclipses 
of (prolate) ellipsoids of different sizes and shapes, but equal in 
luminosity, can be almost exactly represented by that of a pair of 
similar ellipsoids. 

The greatest observable effect will occur when the smaller star is 
much the brighter (Z, small). Even then, the residuals will be very 
small. To express them in terms of the light outside eclipse they 
must be multiplied by 1/7{Ljbi(u, — u.)}. Ina “special Algol sys- 
tem”’ of this type it will be very rare to have b, > 0.4, or u; > 0.3, 
which makes the foregoing expression 0.024 at most. The largest 
residuals, near the beginning of eclipse, amount therefore to o"0o6 
or less. Those close to totality, when expressed in magnitudes, may 
be larger; but it is to be remembered that the adjustment of the 
curves has been made in a summary fashion, and a complete adjust- 
ment, on least-square principles, would give a better fit. 

For k = 0.6 the residuals are still smaller. The residuals indicate 
that the light of the dissimilar pair 1 near the beginning of the eclipse 
will be fainter than that of the similar pair 3. (It is brighter than 


ik 
a 


TABLE 5 


EFFECTS OF DISSIMILAR SHAPE 
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that of the intermediate pair 2, but the final adjustment reverses 
this in this region.) 


COs ¢: a p x y 2 Res. aC, Res. D-—U 
k=0.8 
+1.00 ©.000]+1.000]| 0.00} 0.00] 0.00]} ©.00] 0.00]/ 0.00} 0.00] 0.000 
+0.90 0.033|/+0.821]| 0.17] 0.72] .35 .14/— .25]— .020 
+o.80 0.091|+0.647]| 0.25} 0.91} .50 30I+ .14 .34/— .016 
+0.60 ©.249/+0.310]] 0.41] 1.00] .70 .00 .67 ere) .000 
+0o.40 0.413} ©.000]] 0.59] 0.92] .82 .10 .10]+ .004 
+0.20 ©. 580]—0. 273|] 0.80} 0.7 .go .61|— .06 .07|/+ .002 
0.00 ©.723|—0.500]} 1.07] 0.60] .93 .47 .00 .58 .00o/— 
—0.20 © .834/—0.673]| 1.36] 0.45] .90 .05 .39|/— .05|— .002 
—0.40 ©.909|—0.800]| 1.63} 0.33] .82 .04 .231—  .04 .000 
—o.60 0.965|—0.8go}] 1.78] 0.23 70 .06 .00 
—o.80 0.987/—0.952|| 1.61] 0.14] .50 . 06 .03|-+ .06/— 
—0.90 ©.995|—0.978]| 1.28] 0.10} .35 .o7||-+ .o1j+ .o7/— .002 
—1.00 I .000]—1.000]} 0.00] 0.00] 0.00]/ 0.00} 0.00]/ 0.00} 0.00} 0.000 
k=0.6 

+1.00 0.000]-+1 .000}]} 0.00] 0.00] 0.00]| 0.00] 0.00} 0.00] 0.000 
+o.90 0.029|+0.849 10] .68] .26 .o9]— .12]/— .016 
+o.80 0.081|+0.688 15| .85| .37 .08 .22}— .og|— .o10 
+o.60 212|+0. 395 -231 .50 .45 .00 .000 
+o.40 ©. 361|+0.125 .58 -41/— .04 .04/+ .007 
+0.20 0.508}—o.118]} .63 .40/— .03 .03/+ .004 
0.00 0.644|—0. 333 .80] .64 ere) .002 
—0.20 0.762/—0.518 .63 .24|-+ .02 -30/— .06/— .004 
=0.40 ©. 856|—0.675 «SS! .13|+ .03 
—o.60 0.928)/—0.805]| .41] .50 .00 
—o.8o 0.977/—0.912|| .48] .27] .37 .04|-+ .003 
—0.90 0.992/—0.951 .36] .19] .o2|i+ .o2/-+ .o5|/-+ .002 
—1.00 I .000]—1.000]] 0.00] 0.00] 0.00]/ 0.00] ©0.00]} 0.00] 0.00] 0.000 


This also happens if we try to represent the eclipse of a disk dark- 
ened at the limb by one of a uniform disk. For the first we have 


sin? 6 = A + ByY(k, a). 


If we denote the values for uniform disks by accents, we have also 


sin? 0’ = A’ + By’(k’, a’). 


. 
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With the tables of the ¥-functions;’ it is easy to make the curves 
cut at any three desired points, and so find the difference a’ — a 
for any given value of a. The results are given in the last column of 
Table 5. The curves are made to cross at a = 0.25, 0.70, and 0.95 for 
k =0.8, and at a = 0.20,0.60,and 0.90 for k = 0.6—close enough to the 
intersections in the case of the dissimilar stars to give substantially 
the same type of adjustment. The difference of brightness, D — U, 
is here given in terms of Z, as unit. The course of the curves of 
residuals is strikingly similar (as it has been practically forced to be). 
It appears that the residual effects of the difference of shape are very 
similar to those which would be produced by a slight darkening at 
the limb of amount approximately equal to twelve and eight times 
the unit of the residuals, which, in the extreme case just assumed, 
would amount to 0.3 or 0.2 in the coefficient of darkening. 

When wu, — uw, is positive (as it is in the special Algol type), the 
effect simulates an increased limb darkening. This is opposite to 
Walter’s conclusion; but he did not take into account the effects of 
rectifying the light-curve or those of adjustment of the elements— 
except to compare the effects of shape, without such adjustment, 
with those found by Dugan for limb darkening after adjustment. 

By introducing such a degree of darkening into pair 3, its light- 
curve can be made to agree with pair 1 far within the accuracy of 
present observations (except perhaps when the minimum is very 
deep). To attempt to detect a difference of shape by photometric 
observations in a single wave length would therefore be almost hope- 
less, if we could be sure in advance what was the degree of darken- 
ing; otherwise it would be effectively impossible. 

It does not seem necessary to compute additional cases--espe- 
cially in view of the excellent agreement of these two. 

For nearly equal radii the problem is more difficult. Equations 
(23) and (24) are no longer applicable if the radius of curvature of 
the larger disk, at the end of its major axis, is less than that of the 
smaller—that is, if b,/d? < b./d3. In this case, the last parts of the 
small disk to disappear are two symmetrical areas on each side of 
the axis. If this condition is nearly but not quite satisfied, the con- 
vergence of differential approximations will be bad. 


37 Ap. J., 35, 335, and 36, 245, 1912. 
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When k’ is small, we find (neglecting &’4 in the parentheses) 
x = sin oi {1 — k’ cos + — 3 sin — cos 


y = sin g3{1 + k’ cos — 3k” sin? , 


| 


= k’ sin + gk”? sin? , 


= k’ sin’ + gk’ cos + — 12 sin? . 


If 
S = I = sin’ 
+ k” sin cos + sin} + 2 sin’ 
mC, — 3S = cos sin ¢4(sin? — 
2- 49S =§ ¢; sin 3(sin? — 1%) . 
If 
x 
~ + k’y + kz 


then, neglecting k’4, 
S’ = 2k” cos sin 
whence 
tC, — 83S — mS’ = &k” cos sin ¢3(sin? — — tm). 
We may thus make the computed curve intersect the given curve at 
any two points equidistant from the middle, as well as at the middle 


(cos ¢, = ©) and at the ends. 
We will then have 


rX = mk’ ; 
= 84 (3+ m)k”, 
or, by (51), neglecting small terms, 
cot i’di’ = = — 
1 : 
db,_ _ 3 (1 + ; _ 3 (1 — 
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These have to be multiplied by b](u, — u.). For very small values 
of k’, they would nevertheless become large enough to vitiate the 
method of differential corrections, but this is of no practical impor- 
tance. 

A determination of the degree of darkening from photometric ob- 
servations—even if they covered the whole light-curve and were of 
the exorbitant precision required—would give a value of the dark- 
ening coefficient increased by the effect described above. 

Differences in shape may, however, be determined by simultane- 
ous precise observations in two well-separated wave lengths, as 
Walter has recently shown.3* If the degree of limb darkening and 
the coefficient of the gravity effect do not change with wave length, 
the problem is simple, for the individual photometric ellipticities of 
the stars will be constant, while the observed values are means 
weighted according to L,, L,. Such differences increase the number of 
unknowns and demand additional data. Very accurate observations 
of the principal minimum should determine the difference of the dark- 
ening coefficient for the bright component, irrespective of the shape 
of the companion, for it appears safe to assume that the geometrical 
aspects of the eclipse are independent of the wave length (except for 
stars with extended envelopes, like ¢ Aurigae). The absolute dark- 
ening coefficient for the companion might be determined, in princi- 
ple, from the secondary minimum, for the radius and luminosity of 
the smaller star may be found from a deep principal minimum alone. 
Sufficiently numerous and precise observations of a star like U 
Cephei or U Sagittae might, therefore, solve the whole problem— 
subject to the uncertainty aroused by the existence of the unex- 
plained asymmetry displayed by some of the best-observed light- 
curves. 

IV. SUMMARY AND CONCLUSIONS 

1. In Part I it has been shown that in all short-period binaries 

for which apsidal motion has been detected or has been proved to 


be very slow, the effective internal concentration of density must be 
high. For the eight systems for which tolerably good values (not 


38 Zs. f. Ap., 16, 167, 1938. 


| 


ELLIPTICITY IN ECLIPSING BINARIES 673 


marked with colons in Table 3) can be obtained the mean value of 
the constant K is 0.011. It follows from equation (3) that the super- 
ficial ellipticity of figure of such stars will be, on the average, only 2 
per cent greater than that of the Roche ellipsoid (eq. [13]). 

2. In Part II it is shown that the photometric ellipticity, derived 
in the usual way from the part of the light-curve outside eclipse, de- 
pends on the internal constitution (K), the degree of darkening at 
the limb (x), and the gravitational effect in surface brightness (y) 
(eq. [15]). The effect of internal constitution (if at all similar to that 
of the stars for which it can be determined) is wholly negligible in 
comparison with those of the others. 

3. In Part III the effect of difference in shape of the two compo- 
nents of a close eclipsing pair are discussed, assuming them to be dis- 
similar prolate ellipsoids of different but uniform surface brightness. 
It is found that the light-curve can be represented within the errors 
of ordinary observations by the eclipses of two similar ellipsoids 
having an ellipticity which is the weighted mean of those of the 
actual components (eq. [33]). The small outstanding differences can 
be almost perfectly represented by a slight darkening at the limb. 

It appears, therefore, that there is little or no hope of determining 
precise geometrical figures of the components of an eclipsing pair 
from photometric observations. Greater accuracy can undoubtedly 
be obtained (when the masses are known) by calculating the geo- 
metrical ellipticity by equation (3) (assuming K = 0.01 if the re- 
finement appears worth while) and utilizing the photometric obser- 
vations for the study of darkening at the limb and of the gravity 
effect. 

In the determination of orbital elements for eclipsing pairs the 
light-curve is usually rectified with the photometric ellipticity con- 
stant z derived from the non-eclipse portion, and the same constant 
is introduced into the geometrical equations for the determination of 
the radii and inclination. It would probably be a better approxima- 
tion to use the value of z computed dynamically (with K = 0.01) in 
these equations. The results will not be rigorous in any case, for in 
the discussion of eclipses of ellipsoidal stars the assumption has 
been made, to avoid great complications, that the contours of equal 
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apparent brightness on the disks are ellipses similar to, and con- 
centric with, the limb. An exact calculation, taking into account the 
effects of gravity as well as darkening, would be very complicated. 

This problem has been discussed by Takeda*’ in considerable de- 
tail. His work deserves careful consideration by anyone who pursues 
the question. Very accurate observations—preferably in widely 
separated wave lengths—will probably be required if definite results 
are to be reached. 

PRINCETON UNIVERSITY OBSERVATORY 
June 13, 1939 


39 Mem. College Sci., Kyoto Imperial U., A, 20, 47, 1937. 
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THE DISSIPATION OF PLANETARY FILAMENTS 
LYMAN SPITZER, JR.' 


ABSTRACT 


An investigation is made of the forces acting on a cylindrical filament of stellar mat- 
ter produced by an encounter between two stars. It is shown that under rather general 
assumptions such a filament will expand under its own internal pressure much more rap- 
idly than it will lose energy either by radiation or by turbulent convection. A similar 
conclusion may be deduced for a thin ribbon-like filament produced by a grazing colli- 
sion. It seems probable that a stellar encounter will simply produce an extended gaseous 
nebula around one or more of the stars involved. 


Since Lyttleton’s? important paper there has been a renewed in- 
terest in encounter theories of the origin of the solar system. The 
papers by Luyten’ and by others which have followed have dealt 
primarily with the stellar dynamics of the proposed encounter and 
have been largely concerned with difficulties of energy and momen- 
tum. 

It is also relevant to investigate the conditions under which a fila- 
ment of stellar matter, presumably torn out of both stars by the en- 
counter, can condense into planets. Early work‘ on this phase of the 
problem involved the assumption that the temperatures of the fila- 
ments were of the same order as those of the solar photosphere. As 
Russell’ has stressed, however, this assumption is certainly incorrect, 
since a filament with as much mass as Jupiter must come from layers 
in which the stellar temperature is at least 10° degrees. This fact had 
also been discussed by Jeffreys,° who showed that the planets would 
expand to radii comparable with the distances of their satellites. 
Jeffreys discussed only equilibrium configurations, however, and did 
not investigate either the rate of expansion or the rate of loss of in- 
ternal energy by radiation. 

* This work was done while the author was holding a National Research Fellowship. 

2M.N., 96, 559, 1936. 

3 W. J. Luyten and E. L. Hill, Ap. J., 86, 470, 1937; M.N., 99, 692, 1939; R. A. 
Lyttleton, M.N., 98, 536, 1938. 

4H. Jeffreys, The Earth, p. 27, Cambridge University Press, 1929. 
5 The Solar System and Its Origin, p. 112, New York: Macmillan, 1935. 


6 M.N., 89, 731, 1929. 
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As will be shown hereinafter, the gases within a filament at a tem- 
perature of hundreds of thousands of degrees will accelerate outward 
and within a few hours will reach the velocity of escape. The fila- 
ment will then presumably dissipate. This fate might perhaps be 
avoided by the radiation of energy at a sufficiently rapid rate. The 
optical depth of any filament with appreciable mass, however, will 
be at least several millions, and the time required to radiate half the 
internal energy of the filament would be at least several months if 
there were no expansion. Since no convection currents with veloci- 
ties less than the escape velocity can dispose of so much energy ina 
few hours, the filament must therefore explode.’ 

In section 1 of this paper it is shown that a filament produced by 
an encounter has, in general, a positive energy, i.e., more than enough 
to dissipate the entire filament to infinity. Section 2 presents the 
formulae for the rate of expansion of a cylindrical filament of circu- 
lar cross-section on the assumption that no radiative loss of energy 
occurs. In section 3 the possibility of convection or turbulence is 
considered, and a solution of the equation of radiative transfer is 
given for the case in which the intensity of radiation is a function of 
the time; it is shown that the rate of energy loss is less than a hun- 
dredth of that required to avert the dissipation of the filament. Sec- 
tion 4 extends the analysis to a ribbon-like filament. Section 5 
discusses the relevance of these results to the encounter theory of 
the origin of the solar system. 


1. It is evident that, if a gaseous filament has a positive energy, 
the filament will expand to infinity (provided y exceeds 4/3) unless 


7 This result is reminiscent of A. W. Bickerton’s theory of partial impact (Trans. 
N. Zealand Inst., 11, 125, 1878), in which the appearance of a nova was attributed to 
the explosion of such a filament. The origin of the planets, however, was also attributed 
to this process. Bickerton deserves mention for what is probably the first detailed state- 
ment of the encounter theory of the formation of the solar system. 

The expansion of such a filament has also been discussed by F. Nélke in his thorough 
investigation of encounter theories (Der Entwicklungsgang unseres Planetensystems, 
pp. 188-191, 1930). His analysis is primarily intended, however, to show that the fila- 
ment will leave the gravitational field of the sun, and he considers the motion of a point 
mass in the gravitational field of the two stars, neglecting both the mass and the inter- 
nal pressure of the filament. Nélke’s investigation is thus concerned with the same dy- 
namical difficulties that Lyttleton and Luyten have more recently discussed, and has 
no direct connection with the present analysis. 
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it loses energy at a sufficiently rapid rate. If a filament is to be un- 
stable, then, it is a necessary condition that the ratio of the internal 
positive energy to the absolute value of the negative gravitational 
energy exceed unity. A lower limit for this quantity is easily calcu- 
lated. 

One may assume that the matter in the filament, assumed to be of 
total mass M, was originally in one of the two stars involved in the 
supposed encounter, and that this matter had originally a maximum 
density ps: and a maximum temperature 7,,. Since, in general, the 
temperature within a star decreases outward much less rapidly than 
the density, the mean stellar temperature of the matter in the fila- 
ment may be taken as not less than $7. (In the standard model, for 
instance, the mean temperature for the entire star is 0.585 times the 
central temperature.) 

The total internal energy W of the filament will then be not less 
than 3kM(35T7.1)/2um., ergs, where » is the mean molecular weight, 
m, is the mass of unit atomic weight, and & is the usual Boltzmann 
constant. If y, the ratio of the specific heats, is less than 5/3, W is, 
of course, considerably greater than the foregoing lower limit. 

The self-gravitational energy —® of the matter in the filament 
will have an absolute value not greater than that of a sphere of uni- 
form density ps, with a radius S equal to (3M /47ps)'/3, for which the 
potential energy is —3GM?/5S, where G is the gravitational con- 
stant. If we use the well-known formula (I.C.S. 87-1) for T3/p ina 
gaseous star, assume that radiation pressure is negligible, and use the 
Eddington mass-luminosity relationship (1.C.S. 84-6), a minimum 
value of W/® is given by 


2/ 
W/® > 0.283 (1) 


where M, is the mass of the parent-star. This is greater than unity 
when M is less than 0.15M,. 

We see that W/® is greater than unity for small masses. If M is 
comparable with M,, a separate investigation is required. Let us 
consider the most favorable case, in which a star is split by an en- 
counter into two spheres of equal mass. From the virial theorem it 
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follows that, if radiation pressure is negligible, W/® for the parent- 
star is o.50—almost twice as large as the lower limit found above 
when M is set equal to M,. After the encounter the original kinetic 
energy will be divided equally between the two spheres. The po- 
tential energy of each, however, will be (3)5/3 times the former value 
of —®, and hence W’/® for each sphere will be 0.79. Unless the stel- 
lar matter consists entirely of hydrogen and helium, the additional 
energy of ionization and excitation would suffice to raise W/® well 
above unity. 

In addition, the encounter itself may be expected to increase 
W/® even further. If the passing star is sufficiently massive and 
energetic to split the original star into two equal masses, it will pre- 
sumably give gravitational energy to each component, decreasing ®. 
Physically this means that, if a star is torn in two by tidal action, one 
may expect each component to be enormously elongated and hence 
to have much less gravitational energy. The internal energy W 
could not be decreased by the encounter, since superelastic collisions 
seem impossible. 

In other words, a gaseous star A in equilibrium has a reserve of 
internal energy equal in magnitude to at least one-half its gravita- 
tional energy. A catastrophic encounter with star B will divert en- 
ergy to A from the mutual motion of A and B. This will result in an 
enormous decrease of the gravitational energy of A, reflected partly 
in the disruption of A and partly in the elongation of each component 
of A. If, moreover, the encounter is a grazing one, the thermal en- 
ergy of A will be increased. As a result, the internal energy will ex- 
ceed the value of the negative gravitational energy, and either or 
both components will explode. 

One may conclude, therefore, that any mass of stellar matter torn 
loose from one star by the passage of another will have more than 
enough energy to dissipate itself to infinity. Radiation of energy, 
however, will decrease the value of W; whether or not a filament 
formed by an encounter would actually disrupt completely therefore 
depends on the rates both of expansion and of radiation. 


2. In order to calculate the rate of expansion, certain assumptions 
of uniformity are apparently necessary; these simplifications should 
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not alter the order of magnitude of the results. One may consider 
an idealized cylindrical filament of circular cross-section stretched 
between two stars, one of which is receding with a velocity V relative 
to the other. Let M, R, and z denote the mass, radius, and total 
length of the filament; M’, R’, and z’ will be used to denote the same 
quantities in units of the solar mass and radius. The masses and 
radii of the two stars are left unspecified, as they play no primary 
part in the expansion of the filament. 

The material in the filament will be acted upon by the pressure of 
the gases in the filament and by the gravitational attraction both of 
the filament itself and of the two stars. Only those forces will be con- 
sidered which act in the radial direction in the filament, perpendicu- 
larly to the axis of the filament. The equation of motion in this sim- 
plified case is 


(2) 


where ¢ is the perpendicular distance from the axis of the filament, a 
is the radial acceleration, p the pressure, and ¢ the gravitational po- 
tential. If we assume that all quantities are functions of r and ¢ 
alone, the mass average of the acceleration, which will be denoted 


R op 0g 


The integral of the first term on the right-hand side of (3) will be 
denoted by d,, that of the second, by a,. If now we integrate by 
parts for d, and substitute 7R?zp,, for M, where p,, is the mean den- 


sity, we have 
R 
2 2 
R? pm J Rpm (4) 


by a, becomes 


P,,, the average of P over the volume, will be less than P as defined 
in (4), since the large values of P for small r are weighted much more 
heavily in P than in P,,. We may therefore find a lower limit for 4, 
by setting P equal to P,,. As shown in Figure 1, T does not vary ap- 
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preciably throughout most of the filament, and we may legitimately 
assume that T is constant. With these two assumptions, we have 


where the symbols have the same meanings as before. If the filament 
has been assumed to expand adiabatically, we know that 


(6) 


Pst 


where y is the ratio of the specific heats. Formula (6) is rigorous only 
for the portion of the filament which was initially at the maximum 
temperature 7, and had the initial maximum density p,;. To an ade- 
quate approximation, however, this formula may be applied to the 
mean temperature and density of the gases in the filament. 

To evaluate ad, we appeal to the virial theorem. This theorem 
states that a gaseous system expands when its thermal kinetic energy 
W is greater than @, where —@ is again the gravitational energy. 
Since in the previous section it was shown that W is actually greater 
than ®, one may deduce that 4, is at least twice as great as that part 
of ad, which arises from the gravitational attraction of the filament 
on itself. It is therefore possible to neglect this gravitational attrac- 
tion in the computation of the rate of expansion; the rate so calcu- 
lated will be too great by a factor of less than 2. 

The gravitational attraction which the two stars exert on the 
filament is difficult to discuss in detail but is not likely, in any case, to 
prevent the expansion. Only the radial component of this attraction 
is effective in this connection, and since the force perpendicular to 
the line connecting the two stars varies as the inverse cube of 3, it 
is evident that as soon as 2’ is much greater than unity such a force 
will be negligible. 

We may let ¢, be the time at which this stellar gravitational field 
becomes negligible and the filament begins to expand. Similarly, 
p. and z, will be used to denote values of p,, and gz at the time /,. If 
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t, denotes the value of ¢ for which the radial velocity outward, 2, 
equals the velocity of escape from the filament, v,,, then we have 


(7) 


If we neglect @,, assume that a, is constant through time, and replace 
v,, by its value for a spherical mass M of radius R, we find from (5), 
(6), and (7) that 


(*) ume ( (8a) 


Substituting M/7R?s, for p., we have 
— ty = 1.70 6M sec , (80) 


where R’, 2’, and M’ again denote R/Ro, z/Ro, and M/Mo, re- 
spectively ; tt denotes the value of ¢, derived on the assumption that 
p is constant with time; and y in (80) has been set equal to 5/3, its 
value for a monatomic gas. We assume throughout that uw is unity. 
If 7., is 10° degrees, M’ is 1/500, p. is 1/10, and other quantities are 
of order unity, é{ — f, will equal approximately three hours. Since 
Vo is of the order of 30 km/sec, the change in R during this time will 
be {Ro; the change in z during the same period will, of course, be 
considerably larger. 

More realistically, we may assume that p is proportional to 1/z 
and that z increases uniformly at a rate V, the velocity of recession of 
the one star relative to the other when ¢ is equal to ¢,. The time origin 
is chosen so that z equals Vt. In this case ¢, must be at least as great 
as the value of ¢ for which the filament ceases to acquire additional 
mass, the filament expanding along its axis as the stars recede. With 
this assumption we obtain from (7) 


t= + y) (9) 


to 


where ff — ¢, is given in (8) above. The quantity ¢, is equal to 
2o/V; if V is 10° cm/sec and zg, is 2Ro, ft. equals 1.4 + 103 sec, or about 
half an hour. 
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3. One may consider next the loss of energy from the filament. It 
will be shown that neither convection currents nor radiative transfer 
can be responsible for an appreciable loss of energy from the filament 
in the short space of a few hours. 

Convection currents are easily disposed of. No current of matter 
can reach the surface of the filament from the far interior before the 
time ff unless the velocity of the current is greater than the velocity 
of escape from the filament. More quantitatively, if the condition 
that (¢f — t)v,, be less than R is combined with (7) and with the 
formula for v,, in the spherical case, this gives roughly the condition 
that a, be greater than 2d@,, which, as we have already seen, is ful- 
filled in all relevant cases. Hence, if such currents are responsible for 
any important transfer of heat, they must leave the filament entirely 
and will promote the dissipation of the filament rather than hinder it. 

A number of approximations reduce the problem of radiative 
transfer to a tractable form in the present case. As a first approxima- 
tion we assume that the filament is static. Then the equation of 
radiative transfer, combined with the equation of radiative equilib- 
rium, gives 


(ew) = (10) 

where 
¢= a (11) 


J is the integrated intensity of radiation, a function of position and 
of time; k is the Rosseland mean absorption coefficient; & is the ratio 
of the energy density of radiation to the total energy density within 
the filament; and c and p are the velocity of light and the density of 
matter, respectively. 

In the present case it will be assumed that o and & are constant 
throughout the filament, which is again assumed to be cylindrical 
and of circular cross-section; J is taken to be a function only of 7, 
the distance from the filament axis. With these simplifications, (10) 
becomes 


(12) 


(2) 1 oJ 
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A solution of (12) is given by 


t) = Aer, (13) 


where J,(x) is the Bessel function of the first kind of zero order, and 
A is an arbitrary constant. From the usual boundary condition that 
J equals 2H, or —20/c-dJ/dr, when r equals R, we find that d is 
determined by the equation 


(14) 


where 7 is the optical depth pkR from the surface of the filament to 
the center and where y is the solution of the equation 


_ 3t 
Ji(y) (15) 


When 7 is very large, y is practically equal to any of the zeros of 
Jo(y). 

For the lowest mode of radiative decay y approaches 2.405 as T 
approaches infinity. This gives the mode of longest half-life even 


Fic. 1.—The upper and lower curves represent 7(r) and J(r), respectively, taken 
along the diameter of the filament. The boundary temperature is assumed to be negligi- 
bly small. 


when the variation of J with the other co-ordinates of the filament 
is taken into account. Values of J(r) and of T(r), or J'/4(r), for this 
lowest mode are shown in Figure 1 for a cross-section of the filament. 

One may define ¢, as the time at which the energy density equals 
one-half the value it had at the time /,, where /, is again the time at 
which the filament started to expand; zero subscripts will again be 
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used to denote quantities at the time ¢,. Only the lowest mode need 
be considered, since the others decay more rapidly. If we combine 
(14) with the condition that A(t, — ¢,) equals /n 2, or 0.693, then we 
find 

2.08R?p.k 


cyt (16a) 


SEC , (166) 


where the asterisk again denotes a value derived on the assumption 
that p is constant. Since & cannot be less than o.1g (its value for 
electron scattering) and M’ is at least 10-3, 7 — ¢, is scarcely less than 
a year. 

If the gases in the filament are to condense into planets, they must 
radiate more than half their internal energy before the time ¢,, and 
the ratio of t, — ft, to ¢; — ¢, must be considerably less than unity. 
When p and & are assumed constant through time, this ratio is given 
by (8) and (16), which yield 


7 


Since k should be at least unity, k/£ may be given a minimum value of 
10. We may equate M to one-half the solar mass exterior to the 
radius at which the density and temperature before the encounter 
were ps, and 7,,, and determine these quantities from the usual 
standard model. If z' and R’ are assumed to be less than to and uni- 
ty, respectively, we find that (7 — ¢,)/(ff — t,) is at least 80,000. 

If a uniform longitudinal expansion of the filament is also con- 
sidered, J’, or J/p, must be substituted for J in (10) and an extra 
term, which for small — reduces to (y — 1) J’/&# must be added to 
the right-hand side to take account of the work done per unit time 
per unit mass in the course of the expansion. The absorption co- 
efficient & will also vary in this case but will, in general, increase dur- 
ing an adiabatic expansion, provided that some formula of the form 
xp/T7/? is approximately valid for k. Hence, to find a lower limit for 
t, — t,, we may assume k to be constant. If¢, is defined as the time at 
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which the total internal energy is one-half of what it would have been 
had there been no radiation emitted since the time ¢,, the analysis 
leads approximately to formula (9), with ¢, replacing ¢,, 47 replacing 
ff, and 6 replacing y, where 6 equals 3y — 4. 

To find an extreme minimum value for (t, — t,)/(t, — to)—a 
value which, in fact, is less than any actually possible—we may use 
(g) in both cases, taking an upper limit for ¢; and a lower limit for 
t,. These limits are reached for y = 5/3 and 6 = o (corresponding 
to y = 4/3), respectively. If ¢f and ¢} are both considerably greater 
than /,, this lower limit becomes 


= 38. 8 


_ 3.27° 
£1/292 (18d) 


Even if V is as great as 10° cm/sec, formula (180) still gives a value 
greater than a hundred, provided values of M, p.:, and 7, are again 
determined from the standard polytropic model for the sun, 
while R’ and k/é are taken to be not greater than unity and not less 
than 10, respectively. Thus, we see that in general the idealized fila- 
ment discussed here is unable to radiate an appreciable fraction of its 
energy in a sufficiently short time to avert complete dissipation. 


4. The preceding arguments have been developed for the case of fila- 
ments with circular cross-section. It has been suggested® also that a 
very close encounter or grazing collision might produce a thin ribbon 
or sheet of stellar matter connecting the two stars. Essentially the 
same analysis may be applied in this case as well. 

There is no need, however, to derive another set of formulae ap- 
plicable to this case. It is evident from (18d) that for a filament of 
given mass M and given length z, the least value of (t, — t,)/(t, — to) 
is attained for maximum radius R. It also follows from general prin- 
ciples that the value of this ratio will increase as the initial cross- 
section of the filament is deformed from a circle of radius R to an 
ellipse with semi-major axis R and minor axis 4, keeping the mass of 
the filament constant. In the first place, it is clear from dimensional 
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arguments that (14) must be a general formula in which y is a nu- 
merical constant of the order of unity, dependent on the geometry 
of the particular case envisaged, and in which 2X is the least dimen- 
sion of the radiating system, which in this case is #. (For a slab of 
thickness 2R it may be shown that y is equal to 7/2.) The optical 
depth 7 will remain constant during this deformation, provided & is 
assumed to be constant; d will vary roughly as 1//; and f — ¢, will 
decrease proportionally to h. 

While the escape velocity will not be greatly affected by this trans- 
formation, the acceleration will be considerably changed. The pres- 
sure will increase as p increases; the accompanying increase in T 
will heighten this increase of P. The acceleration is proportional to 
the increase of pressure inward per unit element of mass, and in the 
course of the deformation it will change proportionally to the pres- 
sure itself. Since ff — f, will therefore decrease somewhat more 
rapidly than h, the ratio (t, — ¢,)/(t, — ¢.) will increase slightly as h 
is decreased. Hence, if this ratio is greater than 100 for any relevant 
circular cylindrical filament of diameter 2K, it will be even greater 
than this lower limit for any ribbon of width equal to 2R. 

Convection may also be ruled out in this case. The acceleration, 
as we have just seen, increases roughly as 1/h. Since the escape ve- 
locity is relatively unchanged by the deformation, the time ¢, — ¢, 
during which the convection currents must be effective, if they are 
to prevent the complete dissipation of the filament, decreases accord- 
ingly as h. Since this decrease of ¢; — ¢, balances out the decrease in 
h, the argument of section 3 may again be applied; i.e., if the currents 
convey enough matter to the surface to produce an appreciable 
effect, their velocity must be greater than the velocity of escape from 
the filament. From these considerations and from the conclusions 
reached in section 3, one may infer that a cylindrical filament will be, 
in general, unable to avert expansion and dissipation, regardless of 
what its cross-section may be. 

It is of interest to note the magnitude of these effects. For a fila- 
ment of length Ro, of width 2Ro, of thickness Ro/50, of mass 
Mo/500, and for which p,.?/3/7T., equals 10 °°, the acceleration is so 
great that less than two minutes are sufficient for the gases to acquire 
the escape velocity of roughly 35 km/sec; during this time the thick- 
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ness of the filament increases, on the average, not more than 15 per 
cent. 


5. The foregoing analyses are based on several approximations, 
but there is a considerable margin of safety in the conclusions. Even 
under extreme conditions there is a difference of two orders of magni- 
tude between the short time required to produce macroscopic veloci- 
ties sufficient to disrupt the filament and the long interval necessary 
to radiate half the internal energy. It is not likely that any or all of 
the various approximations made could introduce so large a factor. 

It is difficult, therefore, to see how the filament as a whole can 
avoid dissipation. Although only the average acceleration has been 
discussed, it is unlikely that any part of the filament can remain 
after most of it has been forced off into space. The radial accelera- 
tion near the center of the filament is small, to be sure, but the veloc- 
ity of escape decreases linearly with 7 for a uniform filament. Fur- 
thermore, the central gases will remain at the highest temperature 
for the longest time and will thus have more chance for expansion. 

Nor is it probable that the outer portions of the filament might 
condense into planets as the filament expanded. Any sphere of gas 
in the expanding filament will itself be expanding in a direction per- 
pendicular to the radius vector of the filament with a velocity di- 
rectly proportional to the radius of the sphere. But the velocity of 
escape from any sphere will also vary directly as the radius. Thus, all 
spheres within the filament have the same chance of condensing, in- 
dependently of their radii. The filament as a whole, however, and 
hence a sphere with a radius equal to that of the filament, cannot 
condense after the time /, has been reached. It follows that, as soon 
as the gases of the filament have reached the velocity of escape from 
the filament, no large portion of the filament can condense, indepen- 
dently of the temperature to which it may fall. 

If, then, such a filament would not be likely to condense into 
planets, what would be the fate of the matter so ejected? Such of it 
as did not fall back into the stars or dissipate into free space would 
form an enormously extended atmosphere of some sort around one 
of the two stars directly involved in the encounter or around their 
possible companions. With an assumed mass of 2: 10°° grams and a 
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radius of 30 astronomical units, such an atmosphere would have a 
mean density of 5.4-10~"’ gm/cm and would be comparable with 
the solar chromosphere. Such an atmosphere is reminiscent of the 
Laplace nebular hypothesis, except that in this case there need be no 
lack of angular momentum. The validity of the encounter theory as 
an explanation of the origin of the solar system rests apparently on 
whether or not a non-uniformly rotating atmosphere could condense 
into solid bodies. 


HARVARD COLLEGE OBSERVATORY 
April 27, 1939 
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PROPER MOTIONS IN THE GALACTIC 
CLUSTER NGC 2548 


E. G. EBBIGHAUSEN 


ABSTRACT 

Proper motions are derived on Yerkes plates for an area 30’ in diameter centered on 
the cluster NGC 2548. The results are contained in Table 2. Plots of the motions are 
found in Figure 1, A and B. The positions of the probable cluster members are found 
in Figure 2, A and B. Figure 3 is the spectrum-magnitude diagram of the cluster. The 
luminosity function of the cluster is given in Figure 4. 

This paper gives the results of the determination of proper motions 
in the cluster NGC 2548 (a = 8'08™8, 6 = —5°30’, 1900; 1 = 106°, 
b = +17°). For this cluster Dr. R. J. Trumpler gives a distance 
of 470 parsecs, a diameter of 30’, and the type is designated as 1~2a, 
indicating that a small number of giant stars may belong to the 
cluster. Verification of the membership of these giants is one of the 
main objects of this paper. 

I. THE MATERIAL 


The material available for measurement consisted of four pairs 
of plates, taken with the 4o-inch refractor of the Yerkes Observa- 
tory, each pair having only one image. Table 1 gives a list of these 
plates. 


TABLE 1 
THE MATERIAL 
Limiting 
Old Plate New Plate aii Photovisual Symbol 
in Years 
Magnitude 
F 500 27.9 12:8 PI 
F 498 21.0 11.6 P Il 


II. MEASUREMENT AND REDUCTION 


The procedure adopted in this paper is the same as that used in 
the measurement and reduction of the material for NGC 752." Brief- 


‘Ap. J., 89, 431, 1939. 
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ly, that procedure is as follows: Each old plate was matched by a 
new plate taken through the glass. The two plates were placed film 
to film, and the measurements of Ax and Ay were made in both the 
direct and the reversed directions. The measured values of the dif- 
ferences were then averaged for each star, to form a mean value for 
both Ax and Ay. The measuring device was the same as that used 
for NGC 752. The reduction constants were obtained by a least- 
squares solution of the formula 


Ax =ax+by+c, 


with a similar formula for the y-co-ordinate. 

In order to select the stars of small motion for the least-squares 
solution, a preliminary reduction was first made for P I, the plate 
with the faintest limiting magnitude. Of the 114 stars on this plate, 
76 were selected as comparison stars; for the remainder of the 
pairs—P II, PIII, and PIV-—the comparison stars were always 
taken from this group. The smallest number of comparison stars 
used was 38 for P III. 

Then for each plate and each star the computed values of the dif- 
ferences (Ax, and Ay.) were obtained, and the proper motion in 
both co-ordinates (u, and uw,) resulted from the difference between 
the observed and the computed value of each difference. With the 
aid of the scale value (1 mm = 1076) and the interval in years, the 
values of u, and yw, were converted into units of ten-thousandths of 
a second of arc per year. 

In order to select cluster stars for the elimination of the magni- 
tude equation, the measures of all stars common to all four plates 
were averaged and plotted. The resulting plot indicated a very dis- 
tinct concentration of stars near the origin, and from these stars 
41 were selected as provisional cluster stars. The plate and filter 
combination used for the plates of this cluster give magnitudes very 
nearly on the photovisual system; but since no photovisual magni- 
tudes were available, the measured diameters were used in the re- 
ductions. Subsequently, these diameters were converted into ap- 
proximate photovisual magnitudes by means of the A-type stars 
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and their photographic magnitudes. For only two plates, P I ys 
and P III ws, was any magnitude equation detectable. 

After the magnitude equation had been removed, the material was 
ready for the determination of the weights of the individual plates. 
Again the same plan was followed as in the case of NGC 752. The 
differences PI — PII, PI — PIll, PI — PIV, and Pll — PIV 
were formed for both uw, and w,. Each set of differences for each 
plate was divided into three subgroups, depending on photovisual 
magnitude. These groups were 8™5-9™5, 9™6-11™2, and 11™3~-12™2. 
Then, for each of the twenty-four subgroups the quantity 2v?/n — 1 
was formed. None of these values of the mean error differed signifi- 
cantly from the mean to indicate a real difference; as a conse- 
quence no distinction in weight was made with reference to image 
diameter, and each plate was given unit weight. 

With the aid of these weights the motions from the various plates 
were combined and are given in Table 2. The columns give, re- 
spectively : 

1. The author’s catalogue number (no previous catalogue has 
been published). 

2 and 3. The co-ordinates, on a system in which +x and +y in- 
dicate increasing right ascension and declination, respectively. One 
unit = 876. 

4. The photographic magnitude in Trumpler’s system.? Most of 
the magnitudes were kindly supplied by Dr. R. J. Trumpler in ad- 
vance of publication. For the remainder, interpolations were made 
on an Eastman 4o plate of the cluster taken with the 24-inch re- 
flector of the Yerkes Observatory. For a few faint stars a slight ex- 
trapolation was necessary. 

5. The spectral type, kindly communicated by Trumpler.? A dash 
in place of the spectral type indicates that Trumpler supplied the 
magnitude but gave no spectral type. 

6 and 7. The proper motions in x and y expressed in units of 
0001 / year. 

8. The weight of the motions in columns 6 and 7. Unit weight 
corresponds to a mean error of +0%0024 as obtained from the inter- 


2 Private communication. 
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TABLE 2 
THE CATALOGUE 


No. x y Mog Spec. My Wt 
246 103 + — 85 4 4 
Ace 215 116 + 25|— 81 4 4 
243 117 + 41] + 20 3 I 
225 123 + 9] + 16 4 I 
219 183 — — 18 4 3 
IQI 66 9.0 | gG5 — 1/+ 4 4 I 
199 112 A + o| + 6 4 I 
ere 194 118 E253 — + 4|+ 2 2 2 
205 126 22 | + I0 3 2 
207 132 — 63 | + 22 4 4 
188 14! 9.0] A2 6] + 23 4 I 
“200 164 — 16 | — 10 4 I 
166 46 10.4 Ao +11/+ 6 4 I 
ee 172 48 10.8 Ao + 12] + 12 4 I 
171 60 | = 20 I 4 
Oe 167 90 Poe Go —152 | + 48 4 4 
173 121 12.2 — + 22/+ 7 2 2 
177 133 Wes — 74/+ 3 I 3 
176 144 A — + 12 4 I 
178 139 + 14] — 63 I 3 
166 151 10.0 B8 4 I 
182 205 + 86] + 69 4 4 
152 32 9.9 A2 —- 3 3 I 
ee 154 38 II.4 — — 238 | +129 4 4 
ach 158 44 — 41 | — 49 I 3 
159 55 11.6 +13|/+ 2 2 2 
148 63 10.7 A3 7 4 I 
161 70 |, — 11] — 7! 4 4 
154. 83 A + 4 ° 3 I 
160 89 Az + 7/+ 2 4 I 
162 04 — + 21 | + 19 3 2 
150 98 — — 75 | — 13 3 4 
154 106 11.6 | <A3 + 6 3 I 
157 106 10.6 | A2 + 8 ° 4 I 
142 119 + 5|— 28 3 2 
140 5 AS + 30] + 10 3 2 
138 129 8.8 | Az2 + 13 2 I 
142 133 9.7 | Ar + + 14 4 I 
* No. 44: double star, brighter preceding component. 
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No. x y Mog Spec. My My Wt. 
148 132 11.6 A: — 7 3 
149 144 9.2 Al + 11 4 
159 143 — 66] + 19 I 
161 178 — 56] — 53 I 
139 181 8.7 Ax + 2]+ 10 4 
118 48 52.0 — —144 | + 18 3 
121 50 11.6 | As: + 18] — 34 3 
118 54 11.9 — 21 ° 4 
127 63 TE. 7 A2 — 5§ 3 
116 68 11.4 AS + o| + 4 2 
119 72 10.2 | G8 1] 15 3 
128 77 12.2 A: + 19 | — 45 I 
121 84 12.41 + 2 — 37 4 
129 87 10.9 | At —- 1|/- 17 4 
116 96 12.4 — + 15 | — 85 2 
113 97 12.6 + 10} — 16 I 
120 103 12:7) — 2 + 54 2 
134 102 8.9} Ar + 14]/+ 2 4 
BO 135 107 9.9 | gGo + 4 4 
128 10.9 A3 + 25|-— 5 4 
127 118 F2p | — 14]| — 4 
132 122 + 5|-— 14 3 
$33 122 10.5 A + 14 3 
121 10.1 A2 — 8 4 
134 135 + 60 | — 92 I 
III 138 12.5 — — 34] — 24 I 
127 147 10.0 | A2 + 7 4 
132 152 9.9} Ao — 8 4 
131 192 10.4 s 4 
114 197 — 10] — 69 I 
131 232 — 43 | — 67 I 
98 g2 Ag + 10 | — 12 3 
104 04 9.8} Ar +18/+ 4 4 
107 105 “Al + 10 | + 16 4 
97 125 9.6] At I 4 
99 126 —159 | — 58 I 
99 134 52-0 — + 31]/— 5 3 
103 146 10.9 | Go — 238 | —1I10 4 
87 154 | A; + 29|/— 7 4 
87 174 8.4 | Ar — 11] + 6 4 


Member- 
ship 


t No. 69: according to Trumpler, the F2p spectrum is probably composite A+G. 
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TABLE 2—Continued 


No. y Mog Spec. My Wt. 

07 192 —- 1 I 2 
a 103 201 12.0 — g|— 20 2 2 
106 205 + 78] — 30 I 4 
Q4.. go 207 + 17 4 I 
ost. go 209 8.8 At + 8|/— 13 4 I 
96 220 + + 13 2 2 
97. 104 220 12.0 + . I 2 
98. | 82 46 10.6 \o + 4]-— I9 4 I 
99 | 66 71 9.8 \2 —- 8/-— 4 4 I 
100 | 78 116 10.1 \2 + 12/-—- 9 4 I 
IOI 64 128 + 15 | — 27 4 2 
102 67 12.0 — $0 142 2 4 
103 71 172 Oni Al —- 2/— 6 4 I 
104 68 1 11.9 | + 6 2 2 
| 80 212 11.0 2 4 
106... | 37 66 10.7 - — 36] — 9 4 3 
coy 54 85 A2 + 12|- 4 I 
108... 2 86 + 59 | — 21 4 4 
109... 31 106 11.7 93 | 2 4 
49 127 15:7 + I — 93 4 

| | 

1) 7 157 11.8 + 18] — 17 3 I 
| 50 191 12.4 | + 13] + 50 3 
12.4 | + + g2 I 4 
SEA... | ome 194 10.9 — 86| — 45 4 4 
53 202 9.3 + 10/+ 5 2 2 


t No. 95: double star, north following component. 


nal agreement between the plates. A number of faint stars on P I 
did not occur on any of the other three plates; and, although no di- 
rect determination of their weight could be made, it was considered 
reasonable to give them the same weight as the remainder of the 
stars on P I, i.e., unit weight. 

g. A measure of the probability that the star is a cluster member. 
The figure 1 indicates that the star has a high probability of being a 
cluster member; 2 indicates that the probability is lower but that 
the star has a reasonable chance of being a member; 3 indicates that 
the star is probably not a member; and 4 signifies that the star is 
definitely not a cluster member. This assignment of membership 
will be discussed in the next section. 
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III. DISCUSSION OF THE RESULTS 


The data given in Table 2 are complete down to photovisual mag- 
nitude 12.0, within a radius of 15’ of star No. 72. As a result of the 
close proximity of the image of a fainter star to the image of No. 44, 
it was possible to measure its motion only in the y-direction. As a 
result of its motion (us = 0%0013/year) and its position in the spec- 
trum-magnitude diagram (log JT, = 4.04 and myo) = 8.2), No. 44 
may provisionally be considered as being a cluster member. 


+7005 —Toos —"o10 +’o10 +005 —"o05 —"o10 
T T +8 T T T T +6 T T 
A Pe 3on4 B Pa 
bad 
5 staas 2 staas 
oatside. out side. 
Fic. 1 


The weights assigned to the individual stars range from 1 to 4. A 
weight of 4 corresponds to a mean error of +0%0012/year in either 
co-ordinate. Figure 1, A, is a plot of all stars with weight of 3 or 4, 
and in Figure 1, B, are plotted all stars with weight 1 or 2. 

The determination of the probability of membership in the cluster 
was made in the following manner. The average weight of the stars 
in Figure 1, A, corresponds to a mean error of +0*%0013/year. 
Three circles, of radii 2, 2/2, and 4 times the mean error just given, 
were drawn with their centers at the center of gravity of the con- 
centration. Stars within the first circle were assigned to class 1, stars 
between the first and second circles to class 2, those between the 
second and third circles to class 3, and the stars outside of the third 
circle to class 4. 
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The average weight of the stars in Figure 1, B, corresponds to a 
mean error of +07%0021/year. About the center of gravity two circles 
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were drawn with radii 2 and 4 times this mean error. No class-1 
membership was assigned to any star in this group. Class-2 mem- 
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bership was assigned to the stars inside the first circle, and class-4 
membership to those outside the second circle. 

The number of stars with a membership class of 1 or 2 is 74. 

Figure 2, A, is a plot of the position in the cluster of all stars whose 
weight is either 3 or 4, and Figure 2, B, is a plot of all whose weight 
is either 1 or 2. In both figures the filled circles correspond to stars 
with a class membership of 1 or 2. 

Figure 3 is a plot of the spectrum-magnitude diagram. The filled 
circles represent stars with a membership class of 1 or 2, and the 
open circles are those of class 3 or 4. The abscissas are log T., and 


No. of 
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20 
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Fic. 4.—Luminosity function for the cluster members 


the ordinates are apparent bolometric magnitudes. The reduction 
of Trumpler’s spectra and photographic magnitudes to this new sys- 
tem was carried out with the aid of the table given by Kuiper in his 
discussion of the hydrogen content of clusters. 

Figure 4 gives the luminosity function of the cluster with the 
photovisual magnitude as the abscissa. The diagram is complete for 
the area with a radius of 15’ around star 72 and down to about 
My = 12.0. The distance modulus of NGC 2548 is about 8™6.4 


IV. THE PROPER MOTION OF THE CLUSTER 
Since there are only a few faint stars in the cluster, and because 
none of the measured stars occur in either the Boss General Catalogue 
or Schorr’s Kigenbewegungs Lexikon, no attempt has been made to 
determine the absolute proper motion of the cluster. 


3 Ap. J., 86, 176, 1937. 4 Trumpler, Lick Obs. Bull., No. 420, 157, 1930. 
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The writer wishes to express his sincere thanks to Dr. Trumpler 
for having communicated in advance of publication his spectral 
types and magnitudes for this cluster. Grateful acknowledgment is 
made to Dr. Struve for having enabled the writer to continue this 
investigation after he had left the Yerkes Observatory to teach at 
Wilson College. Also, I wish to express my thanks to Dr. Kuiper 
for his general guidance. 
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THE ULTRAVIOLET SPECTRA OF A AND B STARS* 


OTTO STRUVE 


ABSTRACT 

This paper contains (a) the measured wave lengths of 350 absorption lines in the 
A stars a Cygni, a Lyrae, n Leonis, and at Canum Venaticorum, between A 3227 and 
\ 3957; (6) the measured wave lengths of 198 lines in the O and B stars ro Lacertae, 
r Scorpii, 8 Cephei, y Pegasi, and 55 Cygni. The identifications of the star lines with 
laboratory lines are given for each group of stars. 

The Cassegrain spectrograph of the McDonald Observatory is 
equipped with two quartz prisms made by Bausch and Lomb and is 
primarily intended for work in the near ultraviolet region of the 
spectrum. The specifications of the prisms are as follows: 

First Cornu prism: Angle, 62°11’ 
Length of edge, 80 mm 
Length of face, 130 mm 
Second Cornu prism: Angle, 62°11’ 
Length of edge, 80 mm 
Length of face, 140 mm 


The prisms are set at minimum deviation for \ 3933; the total 
deviation is go° 

The collimator has a focal length of 1 m and consists of two 
pieces of UV glass—UBKsg and UZKs, respectively. Its aperture is 
76 mm. The focal length of the camera lens is 50 cm. It consists 
of four pieces of UV glass, having an aperture of 86 mm. Both 
lenses were designed by Dr. F. E. Ross; they were made by Mr. 
J. W. Fecker at Pittsburgh, Pennsylvania. Two additional lenses, 
of quartz, are not yet available. The ultraviolet limit of the trans- 
parency of the UV lenses is about \ 3250. At A 3933 the linear dis- 
persion with the UV camera is about 40 A/mm. At X 3250 it is 
about 20 A/mm. 

The mechanical parts of the spectrograph were designed by Dr. 
G. W. Moffitt and were built under his direction by Mr. Charles 
Ridell in the instrument shop of the Yerkes Observatory. A descrip- 


* Contributions from the McDonald Observatory, University of Texas, No. 13. 
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tion of the instrument by Dr. Moffitt will appear in No. 1 of the 
Contributions of the McDonald Observatory. 

In May and June, 1939, the spectrograph was used to obtain 
nearly 100 widened spectra of bright A- and B-type stars having 
sharp lines. The observations were made in collaboration with Dr. 
A. Unsold, who will use the material for a study of line intensities. 
The following plates were chosen for measurement: 


a Cygni (A2 supergiant) ...Four plates taken on May 26, June 1, and June 2 


a Lyrae (Ao dwarf) ....... Four plates taken on May 25, June 1, and June 2 
n Leonis (Ao giant)........Two plates taken on May 26 and 27 

a Can. Ven. (Ao dwarf) ...T wo plates taken on May 26 

10 Lacertae (Og dwarf). ...One plate taken on May 26 

7 Scorpii (Br dwarf) ...... Three plates taken on May 24, 25, and 26 

B Cephei (Bo dwarf).......One plate taken on May 27 

y Pegasi (B2 dwarf)....... One plate taken on June 1 


55 Cygni (B2 supergiant)...One plate taken on June 1 


All plates are on Eastman Process emulsion. During measure- 
ment two similar spectra were kept in the machine simultaneously, 
although only one was actually being measured. The other plate 
served as a comparison spectrum and was used to verify the reality 
of very faint lines. 

Tables 1 and 2 contain the results of the measurements. The work 
on the A stars is a continuation of Morgan’s study of the spectra of 
A stars in the ordinary photographic region.’ My measurements 
overlap his near \ 3900, and the agreement is very satisfactory. 
The near ultraviolet region of A and B stars has already been in- 
vestigated by Swings and Désirant,’ who used plates taken with the 
Yerkes autocollimating spectrograph at the Perkins Observatory, 
and by Marshall, who measured spectrograms taken with the light- 
glass spectrograph at Ann Arbor. The spectrum of a Cygni has 
been completely measured in the ultraviolet region by W. H. Wright! 
and A. B. Wyse,’ at the Lick Observatory. The agreement of my re- 
sults with those of Wyse’s compilation is very satisfactory. My 
identifications and assignments of the major contributors to the lines 
were made quite independently, and there occur a few—rather un- 

t Pub. Yerkes Obs., 7, Part III, 1935. 

J 83,35; 1035: * 4 Lick Obs. Bull., 10, 100, 1921. 

3 Pub. Obs. U. of Michigan, 5, No. 12, 1934. 5 Lick Obs. Bull., 18, 129, 1935. 
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important—differences in the assignments of the various laboratory 
lines. My plates show considerably more lines than those listed by 
Wright, as should, of course, be expected in consequence of the fine 
optical quality of the new Cassegrain spectrograph. Wyse has al- 
ready pointed out that a study of the spectrum of a Cygni with rela- 
tively small dispersion cannot be considered as final. But the star is 
so important in many spectrographic investigations that a fairly 
complete list of lines observable with our instruments will be valu- 
able. Moreover, the lines of a Cygni are by no means narrow. They 
look as though they are appreciably broadened by turbulence. 
Hence, the blending of lines is essentially intrinsic in character. 
Relatively little will be gained by the use of larger dispersion.® The 
ultraviolet spectrum of a Lyrae has been investigated by R. W. 
Shaw’ from plates taken with a quartz spectrograph at the Lowell 
Observatory. In the region \ 3000 to 3480 the wave lengths of 118 
lines were determined, but no lines were seen between \ 3480 and the 
limit of the Balmer series. My plates confirm the general weakening 
of the lines in the region of strong continuous hydrogen absorption, 
in a Lyrae and in other A stars. But some 40 lines were measured in 
the interval AX 3480-3647. The equivalent widths of these lines will 
yield quantitative information concerning the amount of continuous 
absorption at different wave lengths. 

The ultraviolet spectra of the B stars have been recently investi- 
gated at the Mount Wilson Observatory: Adams and Dunham? 
have given lists of wave lengths ranging from \ 3023 to A 3664, for 
several stars; Sanford? has measured the ultraviolet and violet 
spectra of + Scorpii and has identified several new Ne 1 lines. The 
new measurements extend my previous work in the photographic 
region of the spectra of B stars'® and partly overlap with the work of 
Kihlborn." 

In making the identifications Miss Moore’s A Multiplet Table of 
Astrophysical Interest (Princeton, 1933) was extensively used. I am 


® An extension of this work with the coudé spectrograph of the McDonald Observa- 
tory is being planned. 

7Ap. J., 82, 87, 1935. 

5 Ap. J., 87, 102, 1938. 10 AD. 225; 1031; 321; 1033. 

9 Pub. A.S.P., 50, 244, 1938. ™ Veréff. Berlin-Babelsberg, 12, Part I, 1938. 
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very deeply indebted to Miss Moore for unpublished lists of im- 
portant lines of Fe u, Cr u, and V u. To Dr. P. Swings Lam grateful 
for a long list of unpublished wave lengths of Fe 1u. The original 
literature was frequently consulted, and many additional identifica- 
tions could be made by going to the sources. 

The intensities of the star lines in Tables 1 and 2 are rough esti- 
mates and are intended solely to serve as a guide in making proper 
identifications. Intensity 1 means that the line was visible on the 
two plates kept in the measuring machine. These lines are regarded 
as certain. Intensity 1-o means a line which was clearly visible on 
only one plate. The vast majority of these lines is real, but a few 
may be spurious. Intensity o-1 indicates a line which is visible on 
one plate and which is regarded as suspicious. The letters “n,” 
“nn,” “nr,” and “nv” mean, respectively: nebulous, very nebulous, 
nebulous and shaded toward the red, and nebulous and shaded 
toward the violet. 

The precision of the measurements varies greatly with the in- 
tensity of the lines. All lines of intensity 2 and stronger are meas- 
ured to within a few hundredths of an angstrom unit in a Cygni, the 
star which yielded the most precise results. Very faint lines, or nebu- 
lous lines, may be in error by a large fraction of an angstrom unit. 
Among the B stars the most trustworthy wave lengths are those 
of 7 Scorpii. But even these may occasionally be in error by as much 
as o.5 A, or even more. 

The wave lengths of the laboratory lines are essentially those of 
Miss Moore’s table. Accordingly, solar values are given where 
these are listed by her, and laboratory values are given only for 
those elements which are weak or absent in the spectrum of the sun. 
The intensities of these lines are those determined in the laboratory. 
They are, of course, comparable only for lines of any given atom or 
ion. The most important contributor or contributors of the star 
lines are designated by the plus sign (+) for a Cygni. Caution 
should be used in applying these results to other A stars. For the 
B stars the contributors vary so much from star to star that no at- 
tempt has been made to assign them. Some laboratory lines which 
appear to be too weak to account for the star line, or which differ 
appreciably in wave length, are given in parentheses. 
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TABLE 1 
A STARS 
a CYGNI a LYRAE n LEONIS a CAN. VEN. IDENTIFICATION 
Con- 
N I N I N I r 1 | Ele r I | tribu- 
ment 
tor 
Ti II 3229.21 40 + 
Ti 3229.43 
Tim | 3232.29 | 30 ss 
3234.52 
Ti 3230.50 70 
Feu 3237.82 8 + 
Tiu 3239.05 60 
Tiu 3239.67 
Ti 3242.00 | 60 + 
3247.21 6.94 Fei 3247.17 9 + 
Feu 3247.39 
Feu 3249.91 
Tiu 3251.90 | 30 + 
3258.87 6 9.31 8.92 Fei 3258.77 10 
Crit 3258.77 
Fe 3259.05 10 + ¥ 
Cri 3264.26 12 + 
Feu 3267.04 
3268.09 8.32 Vir 3267.71 |1000 + 
Cru 3268.48 
Fei 3268.51 3 + 
Cri 3269.77 
Cri 3270.14 20 + 
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TABLE 1—Conlinued 
a CYGNI a LYRAE n LEONIS a Can. VEN. IDENTIFICATION 
Ele- Con- 
I I 1 I ment I | tribu- 
tor 
3271.55 4n ok 1.82 I Vu 3271.12 |1200 
Tiu 3271.67 25 4 
Tiu 3272.10 25 
3270.35 3 share 6.04 I 3275.92 3 
Vu 3276.12 |1500 + 
Feu 3276.61 5 + 
3279.25 I 9.39 | I-on|].... 3278.30 30 + 
Fe Il 3279.05 2 + 
Vi 3279.84 | 300 
3281.21 4 1.22 3 .| Feu 3281.29 7 
Vu 3282.53 | I50 
Crit 3283.04 8 
Crit 3285.96 8 
3287.65 = 7.81 3 7.66 2 8.01 2 Feu 3287.47 ay AP 
Tiu 3287.67 40 + 
Cri 3288.04 5 
3289.39 3 9.46 2 9.24 I Feu 3289.35 7 + 
Vu 3289.39 | Ioo 
3201.74 1.86 3 1.85 2 Cru 3201.75 20 
3293.48 pr) 
3205.57 10 5.67 4 5.57 |] 4n 5.64 5 | Fem | 3295.24 4 
Crit 3295.42 30 Tt 
Feu 3295.81 6 
3298.01 8.11 I 8.08 I Fe Il 3297.89 5 + 
Vu 3298.74 | 130 
3303 .19* 2.47 2 Fe II 3302.86 4 
Feu 3303.47 4 
3307.10 3 7.06 2 6.04 2 3300.95 12 
Cri 3307.04 30 + 
3308.76 2 9.15 I Con 3308.15 5 
Ti 3308.82 8 + 
3312.08 2.11 2 2.26 I Cri 3311.93 20 + 
Cru 3312.18 18 + 
Feu 3312.72 I 
* This line is blended in a Cygni but is resolved in n Leonis. 
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TABLE 1—Continued 
a CYGNI a LYRAE 7 LEONIS a CAN. VEN. IDENTIFICATION 
Ele- 
tor 
3314.03 2 4.16 4-37 2n Feu 3314.00 
Crit 3314.06 6 
Cru 3314.57 
3315.33 2 5.52 I 5.01 Cru 3315.20 
Tim | 3315.33 
3317.99 3 7.92 2 Tiu 3318.03 
(Feu 3318.61 | 
Tiu 3321.71 25 
3322.97 6 3.05 4 3.12 4 2.88 6 Cri 3942.66 | F 
Tim | 3322.93 
Fett 3323.07 8 
Ti 3323.40 be 
3324.20 5 4.20 : 4.45 3 4.25 4 | Crit | 3324.05 20} + 
Cri 3324.10 
Cru 3324.34 12 + 
3326.85 3 6.94 2 6.94 In 6.82 I Tiu 3326.78 Ce Aree 
3328.34 2 8.32 I 8.39 I 8.44 I Cru 3328.35 oo See 
3320.51 5 9.42 2 9.62 I 9.47 3 Tin | 3329.5 gl SOEee 
3332.07 4 1.97 2n 2.24 I 2.33 rn Tiu 3332.11 
3335.22 6 5.17 3 5.20] 3 5.41 3 M |\3335.22/ | 4° 
Crit | 3335.28 
Cri 3335.46 
3336. 26 6.52 I 6.49 2 6.38 3 Cru 3336.32 
3337-51 7.64 I Vu 3337.84 
Tiu 3337.85 2 
3338.53 I 8.37 2 8.63 | I-o 8.84 2 Feu 3338.52 ca) See 
33390.908f | 7nv | 9.64f I 40.00f} 4n | 39.90 4 Crit 3339.80 50 | + 
Si ll 3339.84 
Cri 3339.90 
3349. 333 
1\3340.302/| 35] + 
3340.33 \ 
3341.92 5 1.83 3 I.92 Ti 3341.84 
3342.53 5 2.56 2 2.76 2 | 3n | Cru 3342.51 
3346.73 4 6.98 2 6.81 I 6.89 2 Tiu 3346.7 it See 
3347.89 5 I 7.88 2 7.84 3 Crit 3347 .83 
t Blended in a Cygni with following line; also in » Leonis, but is resolved in a Lyrae. 
t Blended in a Canum Venaticorum with 77 11 3341.84. 
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TABLE 1—Continued 


a CYGNI a LYRAE n LEONIS a Can. VEN. IDENTIFICATION 
I I Ele- 
I I r I | tribu- 
tor 
3349.22 Io 9.29 10 9.35 6 9.19 4n | Tin 3348.91 10?].. 
Titt | 3349.00 75| + 
Cri 3349.34 6 
Tin |/3349-39| 
3349.47) 
3353.15 3n 3.53 2n 3.18 I 3.10 4 ‘Gre 4-3953.20 18 + 
Secu 3353-74 
3355.65 6.18 | I-o 5.46 1-o | Feu 3356.26 2 
3357-41 2 7-15 I-o 7.58 I 51 2 Crit | 3357.40 4 
3358.51 5 8.72 3 8.57 4 8.50 3 Feu 3358.25 Tl Cee 
Crit 3358.50 75 + 
3360. 27 4 0.29 I 0.32 3 0.29 3 Feu 3360.10 3 
Cru 3300.30 100 | 
3361.36 9 1.39 4 1.36 4 1.69 2 Tiu 3301.19 125 + 
Vi 3301.51 60 
Cru 3301.77 20 
3366. 23 I 6.24 Tiu 3366.18 8 
Feu 3306.96 
3368 o1 7 7.87 3 8.03 6 7.79 2n | Cru 3367.42 ° 
Cru 3368.05 150 + 
3368.07 9.32 | I-o 9.35 2n | Cru 3368.73 
Cru 3309.05 12 + 
Tiu 3369.22 
Feu | 3369.35 3 t 
3372.00 o-I 2°32 2.39 2n Cru 3372.13 10° 
Tiu 3372.23 10? J 
3372.77 \ 
3372.73 7 2.83 4 2.79 Tiu 3372.86, | 
3374.11 3 4.98§ | I-o 3.95 I 4.50 tn | Niu | 3373.98 
Tim | 3374.35 8 
3370.39 6.44 | o-1 6.45 2 Cru 3376.27 


§ In a Lyrae only Jiu. 


4 
3379.95 6n 80.25 gn 8 
79.93 4n 7 
gn | Cru 
ri 3379.82 60 + 
Tin |/3380.26\ 
(3380.31) 30 + 
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TABLE 1—Continued 
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a CYGNI a LYRAE n LEONIS a CAN. VEN. IDENTIFICATION 
Ble Con- 
I I I I I | tribu- 
ment t 
or 
3381.23 1.38 | I-o I Feu 3381.00 
3382.71 4 2.60 2 2.78 3 2.60 3 Cri 3382.68 66 Fy ices 
3383.78 6 3.80 5 3.85 5 4.00 I Tiu 3383.77 eta Rae 
6.15 In 5.82 I-o 5.61 1-o |(Fel 3386.45 
3387.87 5 7.99 3 7.79 3 7.89 | Crit | 3387.73 3 
Tiu 3387.85 50 + 
Feu 3388.13 
3391.47 4 1.26 2n 1.39 2 I Feu 3391.30 
Cru 3301.42 60 + 
3392.97 3.18 I 3.19 I Vu 3392.66 
Cru 3393.00 25 + 
3394.37 | Onv | 4.56 3n 4.40] 3n 4-45 2 Cri 3393.86 30 
Crit 3394.32 3° 
Tiu 3304.55 
3395-74 2 5.98 5.87 I Feu 3305.34 
3397.17 I-o (Lu ll 3397.06 
Niu 3307.84 
3402.46 3 2.18 I-o 2.36 I 2.29 I Tiu 3402.42 cS) EAE 
3403.30 6 3.49 3 3.43 4 3.31 3 Cru 3403 32 
3405.13 5.22 €.27 o-1 |(Tiu 3404.96 1) 
3407.22 4 7.14 2 7.28 2 7.30 I Tiu 3407.2 a eee 
Niu 3407.32 8 
3408 83 6 8.89 3 8.83 4 8.83 3 Cru 3408.77 150 
10.22 | o-1 10.33 Tiu 3409.82 4 
13.24 1-o | Fe1 3413.14 Es 
3416.08 4 5.91 2 5.73 2 Feu 3416.02 
3419.33 9.71 I-o (P11 3419.25 
(Feu 3420.18 | 
3421.26 7 1.02 2 1.19 5 1.23 3 Cru 3421.20 ya eae 
3422.80 7 2.80 3 2.82 5 2.7% 3 Crit 3422.74 OG Besides 
4.39 1-o |(Feu 3424.17 
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TABLE 1—Continued 


a CYGNI a LYRAE n LEONIS a CAN. VEN. IDENTIFICATION 
Ele- Con- 
tor 
Cru 3428.94 6 
3430.70 30.36 I-o (Fe Il 3430.15 pr) 
(Crir | 3430.42 2) 
3433.36 8 3.40 3 3-33 4 3.30 3 Cri 3433.30 75 
3436.09 5.90 I 6.38 2 Feu 3430.11 
3438.86 9.09 I-o 9.11 I 8.91 In | Matt | 3438.06 
3440.63 I 0.96 I 0.41 I-o 0.72 1-o | Fer 3440.63 
3442.07 8 2.17 4 2.10 5 2.25 2 Mn | 3441.08 30 + 
Fe | 3442.24 3 
3444.23 4 4.32 3 4.09] 2n | 3.84 I Tim | 3444.33 30 
3440.40 2 6.53 I 6.41 In 6.38 1-o |(Ni1 3446.27 100) 
Cou 3440.40 100 
(K1 3446.72 ) 
3448.84 9.39 I (Feu 3448.43 1) | 
3450.98 1.22 2 Crit | 3450.84 3 
Feu 3451. 23 2 
Feu 3451.32 2 
Niu | 3454.17 5 + 
3456.68 3 6.40 In 6.64 In 6.90 In Tiu 3456.40 20 
Feu 3450.93 | 
Cru 3457.62 25 
3400.30 6 0.31 3 0.30 4 0.38 2 Mn | 3460.33 20 
3461.51 2 1.58 2 1.74 I Ti 3461.50 20 
3404.27 In 4.06 o-1 4.10 I 48S 2 Feu 3463.07 I 
Cru 3404.02 2 
Feu | 3464.50 3 
3465.68 3 5.82 I | Tin 3405.65 3 
Niu 3405.65 
3408.64 8.64 3 8.67 I Feu 3468.68 8 
3471.48] 1.23||| 2 0.77 I-o Fei 3470.24 I 
Nit | 3471.35 + 
Cru 3472.07 } 


|| Blended in a Cygni but resolved in 7 Leonis. 
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TABLE 1—Continued 
a CYGNI a LYRAE n LEONIS a CAN. VEN. IDENTIFICATION 
Con- 
I I I 1 | Ele. I | tribu- 
ment tor 
3474.06 6 3.95 4 4.11 4 4.04 4 
3475-47 2 §.53 2 5.57 I 5.29 3 
3477-16 4 6.07 2 7.05 2 7.17 I-o 
3482.85 6 2.87 2 2.83 | 4 2.53 2 
3484.11 4-54 o-! 4.07 I-o 
3485.98 I-o 5.74 I 6.43 | o-1 5.84 o-I 
3488. 66 5 8.69 2 8.62 3 8.68 2 
3491.04 = 0.90 In 1.05 I 0.85 I 
3493.48 8 3.27 I 3.38 3 3.47 I 
3495.70 8 5.55 2 5.82 3 5.60 4 
3497.46 | 6nv | 7.74 2 7.65 2 7-77 3 
3499.82 2 9.78 I |3500.02 I 00.04 I 
3504.86 5 4.82 3 4.72 2 4.59 2 
3507.74 2n 7.89 o-I 745 In 7.81 2n 
3510.79 4 0°. 86 I 0.57 2 0.58 I-o 
3511.76 1.83 I 2.01 I 
3514.01 5 3.90 In 3.97 3 3-74 I-o 


| 
| 
a 
| 
"el 3513.53 30 | 
Niu 3513.04 8 + eee 
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TABLE 1—Continued 
a CYGNI a LYRAE n LEONIS a CAN. VEN. IDENTIFICATION 
Ele- 
I I I ont I | tribu- 
tor 
3517.13 7.89 | I-on]........ Vu 3517.30 
3520.09 3 20.19 I 20.48 I-o | 19.99 I Vil 3520.02 120 
Ti 3520.20 20 + 
Feu 3522.11 pr : 
Cru 3522.15 
3524.47 In 4.29 Ni 3524.54 200 
Vu 3524.71 200 
3520.52 5.98 I 6.30 I Fet 3520.04 20 
3528.96 I-o 9.12 (Cru 3528.26 1) 
3535.48 3 5.42 In 5.34 I 5.61 2 Tiu 3535-41 40 aa 
Feu 3535.63 2 
3538.00 8.67 1-o | Feu 3537.50 I 
Vu 3538.24 50 
3541.34 I.10 I-o Fet 3541.10 15 
Vu 3541.34 50 
Fet 3554.94 40:| + 
Fei 3555.08 pr 
3556.82 3 6.62 I 6.61 Vir 3556.80 |1500 |...... 
Fet 3558.53 30 
3561.20 0.60 1.18 In Vu 3560.50 go 
Tiu 3561.58 3 
3564.65 2 5.08 I 4.96 |1-on] 4.44 I Feu 3564.52 pr + 
Fet 3505.40 
3566.04 aig 6.18 2 Tiu 3565.07 6 + 
Feu 3506.05 
Feu 3566.16 3 + 
Vu 3506.18 BOG 
3567.68 8.55 1-o | Scu 3567 .70| 20 
3567.74) 
3570.01 2 0.03 3 0.13 I 0.06 2 Fet 3570.14 100 
c 8) 
572.3 5 Seu |(3572-4 
3572.33 4 2.55 I 2.59 | o-1 1\3572.s7) | 5° 
3576.67 7 6.30 I 6.57 2 Scu {3576.33\ 
Niu 3576.76 + 
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TABLE 1—Continued 


711 


a CYGNI a LYRAE n LEONIS a CAN. VEN. IDENTIFICATION 
Ele- 
I I I I r I | tribu- 
ment 
tor 
3578.54 I-on 8.58 I Cri 3578.69 200 
Tiu 3578.71 ° 
3581.08 5 1.24 3 1.37 In I.40 In Seu 3580.93 2 
Fei 3581.2 250 + 
3585.30 10 5.27 3 5.35 5 5.53 6 Crit 3585.3 60 + 
Crit 585.5 Ke) 
3587.18 2 7.14 I-o 7.74 I 7.31 I Tiu 3587.15 12 
3589.72 3 9g 89 I 9.55 2 Vu 3589.74 |1000 
3592.12 2 Vu 3592.01 800 & 
3593.62 I I Vu 593.32 600 
Cri 3.50 160 
3505.95 3 6.07 I 5.79 I Tiu 3596.06 Cd ee 
3600.77 o-I Yu 3600.74 300 
3603.70 6 3.74 I 3.62 3 3.85 5 Crit 3603 .64 i Gee 
Crit 3603.80 30 
Cri 3603.86 10 
3605.93 In $.72 I 5.47 I Cr 3605.34 7d 
Fe 3605. 48 5 + 
3608.72 3n 8.73 I 8.91 2 8.97 I Crit 3608 . 66 Ce Rare 
3008 . 87 
3610.60 I 0.40 I-o Yu 3611.05 200 
3013.44 4 3.56 In I 3613.21 20 
Tiu 3613.33 pr 
613.81 
S 
cll 3613.88 100 + 
3614.84 2 4-76 I 4.96 I Feu 3614.87 5 
3618.80 3n 8.69 2 Fet 3618.78 125 
Vu 3618.92 200 
3621.26 5 1.18 2 1.44 2 Vu 3621.20 150 a 
Feu 3621.27 6 + 
3624.84 4 5.08 I 4.56 2 5.02 2 Tiu 3624.84 70 + 
3627.14 Feu 3627.17 I 
3628.68 I 8.11 o-I 8.22 1-o | Yu 3628.71 100 
3631.48 8 I.41 2 1.32 3 1.74 4 Fel 3631.47 | 125 : 
Crit 3631.49 50 + 
Cru 3631.72 25 
3634.98 4.39 I-o Crit 3634.04 5 
Feu 3634.89 5 
3641.31 4 I 0.92 I 1.50 1-o | Fel 3641.22 pr 
Ti 3641.34 100 + 
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TABLE 1—Continued 


| 
a CYGNI a LYRAE n LEONIS a CAN. VEN. IDENTIFICATION 
Ele- 
tor 
3642.78) | 
Cru 3643.22 8 
3645 . 29 
S 
1\3645.34 
3647.76 2 I Crit 3647.40 5 | 
Fei 3647.85 100 | 
3050.40 I 0.53 I Cru 3650. 37 20 
3657.47 8.51 I Cru 3658.19 12 
3659.64 2 9.38 o-1 |. .| 60.50 o-I Tiu 3659.76 60 
Tiu 3662.2 
3062. 26 + 
3664.68 2 5.217 I Has 3664.68 + 
Cri 3604.95 20 
3065.98 2 6.29 3666.10 
Has 3009.47 
3677.84 4 7.62 2 8.03 3 Sn a ee Cru 3077.69 20 + 
Crit 3077.86 30 + 
Cru 3077.93 10 
3079.43 6 9.57 I 9.46 6 ; 3679.36 
3082.78 8 2.97 3n 2.85 8 Hao 3682.81 
3685.22 5 5 5.28 5 5.05 2n Ti u 3085.20 250 
3686.87 9 6.83 4n 6.88 10 7.06 3n Hy 3686.83 
3091.59 10 1.42 5n 1.59 15 2.15 3n 3691.56 
3093.97 o-I 3.85 I Fet 3004.03 20 
(You 3694.20 200) 
3697.18 12 6.97 6 7.28 16 7.36 4n Wh; 3697.15 
3703.78 15 3-74 5n 3-93 | 18 3.93 sn | Hie 3703.86 
3700.07 6 5.90 6.38 Cau 3700.04 IO 
Tiu 3700.22 20 


© Probably Cru in a Canum Venaticorum. 


= 
cit 
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TABLE 1—Continued 
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a CYGNI a LYRAE n LEONIS a CAN. VEN. IDENTIFICATION 
Ele- Con- 
I I I ment I | tribu- 
tor 
Vir | 3709.34 40 
Yu 3710. 29 500 
3711.97 15 1.99 7n 2.07 | 20n 2.30**| 7n | Ais 
3712.88 5 Cru 3712.97 40 
5.32 6 5.38 3 Cru 5.19 20 
Crit 3715.45 18 
Vu 4356.46 19200 1... 
8.16 I Vu 3718.16 
3719.99 3 19.82 2 20.31 3 Fet 3719.95 | 250] + 
Feu 3720.17 
3722.01 15 1.81 8n 2.03 | 20n 1.65 8n Ais 3721.94 
3723.69 o-!I Crit 3723.40 7 + : 
Tiu 3723.01 tr 
Fei 3723.91 pr 
3725.40 2 5.03 I Feu 3725.30 3 
3727.33 4 7.64 3n 7.14 7.35 3 Vu 3727.35 |1000 
Cru 3727.37 10 + : 
Fei 3727.04 50 
3732.93 I Vu 3732.76 | 800 
3733-33 40 
3734-35 15 4.47 | 15n 4.48 | 25n | 4.52 | ron | As 3734-37 + 
Fe 3734.88 
3736.96 6 6.95 6 7-14 7 Cal 3730.92 ao 
Fet 3737-14 15 
3738.42 3 9.03 I 8.42 2 Cri 3738.38 
3741.69 6 1.48 2 1.67 5 1.82 3 Ti 3741.65 
3743.43 2 3.57 I 3.31 2 3.66 3 Fet 3743.37  Soeeine 
Vu 3743.01 
3745.63 4 5.69 3 5.94 3 5.92 3 Feu 3745.36 pr |...... 
Fel 3745.58 
Vi 3745.81 800 
3748.53 6 8.68 3 (Tiu 3748.00 | eae 
3748.27 
(Cru 3748.68 | 
3750.17 15 0.23 20n 0.12 | 20n 0.20 12n Hu 2780.35 + 
Vii 3750 
3755-55 5.23 2 Feu 3755.56 
3757-79 5 8.28 2 8.00 2 8.04 I Tim | 3757.69 30} + 
Fet 3758.25 


** Perhaps blended with Cr 11 3712.97 in a Canum Venaticorum. 


| 
{ 
| | | | 
| 
| 3 
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TABLE 1—Continued 


a CYGNI a LYRAE n LEONIS a CAN. VEN. IDENTIFICATION 
Ele- Con- 
tor 
3759-33 6 9.31 3 9.52 | 7 9.48 3 | Tim | 3759.30 | 200] + | 
Feu 3759.40 6 
3761.45 | 6 | 1.43 | 4 1.44] 5 | 1.77 | 3 | Tim | 3761.32 | 200] + 
Crit 3761.69 4 
Tiu 3761.88 15 
Cri 3701.90 4 
3763.98 4 3.5ott] 2 4.01 4 3.91 3 | Fer 3763.80 | 100 
Feu 3764.11 pr 
3765.55 OF, 5.69 I 98 I Cru 3765.62 8 
Fel 3707.21 80 + 
3779.77 7 0.57 20n 0.77 | 452 0.68 isn | 3770.63 
Vu 3770.97 400 ae 
Tiu 3774-05 | 
3779-53 9.08 | In 9.05 2n | Feu 3779-57 pr 
3781.53 tac I 1.97 I Feu 3781.51 I 
3783.42 3.53 I Feit | 3783.35 4 
3788.11 I 7.56 a Aen 8.35 I Vir 3787.24 150 
Fet 3787.89 50 
3795.01 60 + 
3797-97 20 7.89 20n 8.07 | 30n 7.91 2on | Hie 3797.90 
3806.63 6.46 In 6.88 I Fet 3806.70 fo) 
3813.40 I 3-17 Fet 3812.96 40 
Tiu 3813.40 
Tiu 3814.60 
3815.84 3 5.97 2 5.81 I 6.27 I Vu 3815.38 iol SE 
Fe 1 3815 .85 100 + 
3820.35 4 0.33 3 19.91 | 4n | 20.41 I Fet 3820.44 250 
3821.85 1.76 I 2.10 I Fei 3821.97 pr 


tt Mostly Fe1in a Lyrae. 


i 
r 
| 
| | 
| 
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TABLE 1—Continued 
a CYGNI a LYRAE n LEONIS a CAN. VEN. IDENTIFICATION 
Ele- Cow 
3824.86 4 4.45 2 4.90 3 6. 37 2 Feu 3824.91 4 
3825.75 3 6.11 3 Fet 3825.89 
3827.43 3 7.69 2 7.23 2 7.56 2 Fem | 3827.08 4 
Fei 3827.69 
3829.45 3 9.33 3 9.34 2 Mgt | 3820.37 40 
3832.40 4 2.10 4 2.21 3 Mgt 3832.31 SD 
3835.37 20 5-35 20n 5.38 | 20n 5.27 2on | Hy 3835.39 
3838. 27 5 8.20 4 8.02 4 Feu 3838.04 | eee 
Mgt 3838.30 100 + 
3840.94 3 1.00 2n 0.63 2 35 I Fet 3840.45 80 ao 
Fet 3841.06 80 + 
Feu 3841.36 pr 
3843.24 2.94 I-o Mn tt | 3842.98 I 
3847.77 I 8.15 I Vu 3847 .32 100 
Meu | 3848.24 7 
3849.73 3 9-74 In 9.85 3 50.24 rn | Fet 3849.98 40 
3851.46 o-I |...... (Fel 3850. 83 12) 
3853.74 3 3.11 In 3.86 4 3.65 3 Siu 3853.67 3 
3856.05 5 6.23 $ 6.09 6 6.10 6 Siu 3856.03 8 a 
Fet 3856.38 50 
3860.00 3 59.85 3 59.78 I Fe1 3850.92 | 300 | + 
Feu 3860.11 pr 
3862.61 6 2.70 3 2.62 5 2.60 3 Siu 3862.60 Oot ayes 
3864.03 o-I Feu 3863.41 I 
Vu 3863.81 60 
Feu 3863.94 I 
3865.44 2 5.80 2 5.51 2 5.52 2 Fet 3865.54 on eee 
Crit 3865.50 20 + 
3872.08 4 2.38 2 2.66 | 2n 2.59 I Fe 3872.51 60 
3878.56 jn 8.43 2 8.65 In 9.18 2 Fei 3878.58 
Vil 3878.72 300 
3882.48 2 2.31 In (Mn 3883.28 
3886.61 I Fet 3886.30 
3889 06 20 9.30 25n 8 go | 20n Q.12 in | Hs 


48 
| 
| 

/ 

| | | | | | 
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TABLE 1—Continued 


a CYGNI a LYRAE n LEONIS a CAN. VEN. IDENTIFICATION 
I ‘ I Ele- | Con 
tor 
Vir 3896.16 60 | 
| 8.46 I-o Vu 3809.14 200 | 
3900.56 5 0.70 3 0.51 3 0.50 I Tit | 3900.54 70 | 
3903.36 2 $32 I 4.39 In 3.46 I Vu 3903.27 250 | 
| 
3905.80 4 5.78 2 5.74 3 6.05 3 Crit 3905. 64 18 | 
Feu 3900.04 5 | 
3908. 36 I-o 7.88 I 8.390 I | 
3910.99 | 10.03 I-o 10.45 | 
ae | 
3913.50 5 3.52 3 3.54 Tiu 3913.47 60 | 
3914.26 3 1.07 | 3014.33 250 
Feu 3914.48 2] 
3916.36 2 6.64 | o-1 3916.42 200 | 
3018. 43 I 9.11 I-o 8.54 I (Mn | 3918.32 3) | 
3920.32 I 0.70 I 0.74 I Fet 3920.27 20 | 
? 
3922.88 I 3.16 I-o ace Fet 3922.92 25 
3926.38 I 6.48 In ....]| Mnu | 3926.46 4 | 
| 
| 
3927.95 I Fet 3927.94 30 | 
3930.21 3 | , 0.25 I 0.24 2 Fet 3930. 31 25 
3932.15 2 | Tiu 3032.02 
3033.04 Cau 3933.68 400 | 
| 
| 
3936.01 3 Feu 3935.94 6 | 
3938.58 5nr Feu 3938.29 2 
Feu | 3038.07 
3045.04 2n Feu 39045. 23 pr 
3047.92 I O1 3947 33 10 
Gt 3047.51 7 
O1 3947.61 4 
3052.03 | Vu 3951.07 500 
| | | 
3056.74 2 | Fei 39560. 69 12 
| | 


| 
_ 
| 
— 
| 


Tiu 


Feit 


Cru 


Mgt 


Siu 


Manu 


The limiting laboratory intensity recorded in a Cygni is about 2. The following table gives an ap- 
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NOTES TO TABLE 1 


proximate conversion from laboratory intensities to stellar intensities. 


Lab. 


2 


20 


w 


The limiting laboratory intensity is 2 or 3. 


Lab. 


The limiting laboratory intensity is about 4. 


Lab. 


The limiting laboratory intensity is about go. 


Lab. 
go 
100 
200 


Lab. 
40 
80 


The limiting laboratory intensity is about 2. 


Lab. 


8 


Lab. 
3° 
45 

100 


NO + W 


° 
3 


5 


Very improbable. Only \ 3419.25 may be present; there are no other coincidences. 
Very improbable because JP = 4.33 v. It is possible that \ 3446.72 is present. 


Lab. 
20 
50 


Lab. 
200 


Lab. 
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* 
6 2 4°5 
> o-! 8. 
I 2-3 
* Lab. * 
3 ° 20 2-3 
5 I 50 . 5 
6 I-2 75 5 
AE 
* Lab. * 
I 500.. 2 
2 800 3 
2.. . 
f 
3 
Pu 
Kt 
* 
Cri 
| * 
| 


Cou 
Nil 


Niu 


Lau 
You 
Luu 


OTTO STRUVE 


The limiting laboratory intensity is about 10. 


Lab. * Lab. * 


Lab. * 
200.....- I 


Lab. 
I-2 
8.. 1 

The limiting laboratory intensity is about 100. 

Lab. * 

200....... . 

300.. I 

500... 2 


Uncertain. Only one line present. 
Absent in this region. 

One tentative coincidence. Improbable. 
Improbable. Lab. 200; *1 


There remain very few unidentified lines. The strongest is \ 3845.20, intensity 3 in a Cygni. 


718 
| 
| 
eae é. || The strongest laboratory line \ 3501.73 is probably present. ‘The ionization potential is 17.1 volts. 
Be 
The limiting laboratory intensity is about 3. 
Yu 
ag 
| 
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TABLE 2 
B STARS 
10 LACERTAE 7 SCORPII | B CEPHEI | y PEGASI 55 CYGNI IDENTIFICATION 
| ry | 

3263.40 | (Neu | 3263.44 3) 
3265.43 3 | | Om | 3265.46 | 10 
(Fe ut | 3266.88 | 20) 
96.51 | I-o Het 3296.76 | 1 
3312.62 Ont 3312.30 5 

| | S | 3324-85 | 6 
| Net | 3329.18 4 
| | Our | ae 4 
3334.92 | 3 | 4-75 4 4 | Net | 3334.87 | 10 

3344.66 In 4.27 2n 4 Nell | 3344.44 | 5 
3347.66 I | 7.81 I-o (Fein | 3347.70 | 8) 
| | | | Our | 3350.09 | 4 
3354.70 | 2nn | 4.73 | 2n | 4.38) 2 4.39) 2 4.82) 3 | 3354.52 | 1 


719 

| 

| 

| 
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720 OTTO STRUVE 
TABLE 2—Continued 
10 LACERTAE 7 SCORPII 8 CEPHEI vy PEGASI 55 CYGNI IDENTIFICATION 
3360.68 | | 0.55| 2 | 0.611 1-0 | 61.04) oF Mem | 3360.63 5 
| 3.21 | (Ou 3362.38 1) 
3367.42 | 3 2 Nett | 3367.25 | 6 
Nim | 3367.36 7 
| In (Str | 3370.37 4 
3372.52 | Ca ut | 3372.68 8 
Tim | 3383.76* 
| (O m1 3383.85 2) 
(Out 3384.95 4) 
8.37 
90.11 | I-o go.1I0 2 Ou 3390. 25 8 
3407.13 | I 07.31 On 3407 38 7 
Nu 3408.14 3 
413.70) I-O |....... ds 
3417.69 I 7.45 In 7.95 
3440.10 Ou 3440.39 
3444.13 2 4.31 Om | 3444.10 5 
3447-53 | 2 7.55 3 7-59 3 7-59 3 7-37 3 | Her | 3447.59| 2 
(Bir | 3451.25 | 10) 


3454 
Our 3455.12 5 


* Interstellar. 
t The line in y Pegasi cannot be O 11. 


| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
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10 LACERTAE 
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TABLE 2—Continued 
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SCORPII 


| B CEPHEI 


vy PEGASI 


55 CYGNI 


IDENTIFICATION 


3539.45 


3542.91 


3554.10 


3568.49 


6.04 


10.97 


8.48 


o-! 


I 
81.94 I-0 
7-93 I 
7.38 2 
99.22 I-o 


11.82 Inn 


0.07 Inn 


63.58 I-on 


46 


40 


2n 


66 


09 


I 
75¢ «nn 
84 Inn 
61 Inn 
48 nn 
oo |I-onn 
25 
48 In 

} 
40 I 
32 | 2nn 
15 | 
34 3 
35 o-! 
60 4 


Ele- 
ment 

(Het 3456.7 
Het 3405.89 | 
Ou 3470.42 
OU 3470.81 
Het 3471.78 
Het 3478.95 


| 3486.93 
Het 3487.72 
Su 3407 . 34 
Het 3498. 63 
(Fi 3501.40 
(Fi 3502.95 
(Fu 3503.10 
(Fu 3505.61 
Het 3512.50 
Het 3530.50 


Cal | 3537-75 


Neu 3542.89 
Het ([3554.44\ 
3554.57) 


6) 


t The He line is present in 55 Cygni but is not present in » Pegasi, because of Stark effect. 


| 
} 
| 3 
I I | I I 
78 1 
| 87.09 I 7.48 Inn 7 | 
I 
97.23 3 7-42 In 7 5 
| 
12) 
| 
Inn 30.5% Inn 0.37 2n ° 
| 
| Ge 
I-on r.72 | Inn 2.30 | mm |. 7 
Inn 3.93 | Inn |.. 4.29 | 3n 4 
| | 
| 
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TABLE 2—Continued 
| 
10 LACERTAE 7 SCORPII CEPHEI vy PEGASI 55 Cyont | IDENTIFICATION 
3587.04 | 2nn| 7.05 | 3nn 7.33] 3 | 7-42] (Cum | 3585.83) 2 
| | | Her 42)| 
Ci 3587.68 2 
90.30 2 0.56) 2n 0.31 1 | 0.51 3589.67 4 
| Sim | 3590.46 | 8 
| | | Cu | 3590.87 | 2 
94.10] In 4.36) 1 3593.60 3 
3601 I | v4 | 2 | 01.36 1 (Fei | 3600 93 | 10) | 
| Alt | 3601.62 | 6 
3609.17 In 9.68 I 9.74 I Cul | 3608 96 4 | 
| Cur | 3609.40 5 
12.45| O-I | 32 29 I 12.57 I Alul | 3612.35 4 
3613.54 2 3-63] 4 3.60) 3 3.68 | 4 g-74| 3 | Met | 3613564) 3 | 
| 
3634.27 4n 4.36 4n | 4.16] 7 4.390 | 6 | | Hex 2 
52.22 I | 2.12] 1 2.13 r | 2.07 2 Het I 
| 74.13 | 1S | Has 3673.76 


—_ 
| 
| 
Bk 
| 
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TABLE 2—Continued 


723 


y PEGASI | 55 CyYGNI 


IDENTIFICATION 


| 
| 7 SCORPII | B CEPHEI! | 
| | 
| 
86.10 | rnn | 86.60} 4n | Hx 3686.83 |... 
| 
gl In 92.03 | Inn | ox 64} 5n 
| 
3604.33 2 4.19 4 4 | | Nett | 3604.22 | 10 
3696.76 I 6.77 I 7.52) In 6.67 | 2nn | 6.97 5 (Om | 3695.37 | 5) 
Ay 3697-15 |..-. 
3703.65 3.40 | 3-91 4 | 3703.37 5) 
Hw 3703.86 ‘te 
'3705.00\| 
3705.30 5 4.93 | 5.01] Ion 4.99 | 10 4.98| 4 | Het 3 
3709.75 2 9.75 3 | New | 3709.64 7 
12.04 | 4nn | 11.97 6 ss 3711.97 
3712.88 3 3.40 4 cbs OU 3712.75 7 
| | New | 3713.09 | 10 
3718.12 I 7.91 I Sul 3787-97 7 
3722.07 3nn 1.68 | 2nn 1.87) 4n 1.95 | 5nn 1.87 10 | Hy 3721.94 
| 
3727.35 2 aa 4 7-21) 4 7.55 I 7.19 2 New | 3727.08 9 
| Ou 3727-33 8 
| | | 3732.99) 
| | 
3734.27 4.09, Ion 4.12 |1onn| 4.22] I5 His 3734-37 
| 
| 
3748.08 |1I-on | | 
| | 
3750.25 4nn 49.92} 5n 49.88) 12n | 50.12 | 
|| 
3754.71 4 4.69 
3757.26 3 7.34 I 
3762.36 2 2.45 I 2.52 ee Se 


= i 
| 
ele 
4 
| d 
| a 
= 
) 
) 
| 
| | 
| 
aa 
By 
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TABLE 2—Continued 


10 LACERTAE 7 SCORPII 6 CEPHEI y PEGASI 55 CYGNI IDENTIFICATION 
EI | 
| 
3770.65 5 nn 0.55 | Ion 0.51) 15n 0.65 |15nn| 0.67 18 |(Het 3768.80 | 1) 
(Het | 3770.57 | 1) 
Hu 3770.63 | 
Oul 3774.00 6 
3776.98 I 7.30 2 6.66 | 1-o | New | 3777.16 8 
(OU 3777.00 1) 
85.26 | 1-0 | 84.85 3 Het 3784.88 | 1 
3791.41 3 1.45 2 1.87 I Om 3791.26 6 
| 3791.41 3 
3797 -93 7nn 8.05 | 1onn 7.90) 20n 7.90/15nn| 8.o1 18 | Hw 3797.90 
3806.55 4 06.65) 6 06.60 3 06.46 5 |(Het 3805.75 1) 
St ul | 3806.56 5 
3819.61 
3819.77 9.78) 7n 9.75] 10 9.74] 15 9.65 | 12 | Het 4 
3835.45 Ion 5.82 | 15 nn 5.42] 20n 5.35 5.41 18 |(Het 3833.56 1) 
H, 3835.30 
Sul 3838.32 8 
Vu 3838.30] 5 
51.01 | I-o | 50.56 3850.99 | 25 
(OU 3851.04 | 3) 
(Clu 3851.37 | 30) 
56.47, 1 55-94 I 56.03 3 \(Nu 3855.08 | 2) 
Vu 3856.07 3 
Ou 3856.16 5 
60.39 | o-r | Sim 3860.64 4 
(Clu 3860.83 | 35) 
(Clu 3860.98 | 15) 
62.44 | I-o | 62.81 1 |(Siu 3862.60 6) 
[3867.46 
3867.42 | 2 7.76| 3 7-59] 3 7.66] 3 | 7.50] 4 | Her 2 
{eee 71.67] rn 71.72| 2n 72.07 | 2n | 71.68 4 Het 3871.80 2 
76.39 2 6.20! 3 6.34 3876.19 | 4 | 
Ci 3876.41 2 | 


ge 
| 
AM 
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TALLE 2—Continued 


10 LACERTAE SCORPII B CEPHEI PEGASI 55 CYGNI 

| | 
2.33; 2.31 I 2.97 | 1-0 | Ou 3882.19 
3888. 89 9.12 15 nn 8.79 20n 8.95 | 20n | 9.16, 20 | 3888.64 10 
Hs 3889.05 
12.21} 2 12:27] 3 11.81 I 15,72 | tm On 3911.95 | 10 
3019.87§ | In 9:37 | 2 9.29, 3 8.93 | 3 9.14| 3 | Cu | 3018.08 | 6 
Nu | 3919.00} 6 
Ou | 3919. 28 | 6 
20 37, I 20.77 3 20.64 3 Cu 3920.68 | 8 
24.87. I-on 4.89 2 4.88 2 4.42 I Simi | 3924.44 | 4 
3920.74 Inn 7.31 | 2n 6.93, 3n 6.94 | 3n 6.67 4 Het | 3926.53 I 
(Nut | 3934-41 3) 
Nu 3955 85 | 6 
| 3361.59 8 


§ In 10 Lacertae blended with C 11 3920.68. 
|| Probably interstellar. 


© Interstellar. 
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726 OTTO STRUVE 


NOTES TO TABLE 2 | 


Het Many new Het lines have been found in this region. The members of he two diffuse series extend 
very far in the supergiant 55 Cygni. The highest members become appreciably broad and dif- 
fuse—indicating the presence of Stark effect even in these stars. No Stark broadening is recorded 
in P Cygni (Struve and Roach, Ap. J., 91, 727, 1939). 

Bu Only one very weak line in 10 Lacertae. Improbable. 

Cu Present. 

Cin Definitely present. 

Nu _ Present. 

Nut Present. 

Ol Present. 

Ot Strong in 10 Lacertae. 

Fu Probably present. 

Fut One line probably present. 

Neu Strong and numerous. 

Mgu_ Absent in this region. 

Siu Probably absent in this region. | 

Present. 

Sitv Present. 

Put Absent in this region. 

Su Absent in this region. 

Sut Present. 

Fem <A few weak lines present. Very strong in P Cygni. 

Clit Uncertain. \ 3860.8 may be blended with S 11 3860.64 and \ 3851 may be blended with O 1 
3851.04. These are the strongest laboratory lines. 

Arit The strongest line in this region, \ 3545.84, is absent. 

Cait Suggested by coincidences of two strongest laboratory lines \ 3372.68 (8) and A 3537.75 (7) in 10 
Lacertae. 

Several fairly strong star lines remain unidentified. 
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THE ULTRAVIOLET SPECTRA OF 17 LEPORIS 
AND P CYGNI* 


O. STRUVE AND F. E. ROACH 


ABSTRACT 


The ultraviolet spectrum of 17 Leporis is compared with that of a Cygni. 77 11, 
Secu, and Mnit are greatly enhanced in 17 Leporis. Sit is abnormally weak. 
From the relative intensities of three Mgt lines it is found that the dilution factor 
W > 1074. Other estimates suggest 0.1 > W > 0.01. The observations are, in gen- 
eral, consistent with the theory of cycles. 

One hundred and fifteen lines were measured in the spectrum of P Cygni between 
\ 3265 and A 4026. The Balmer series is followed to n = 23. The triplet diffuse series 
of He 1 is found to extend to appreciably higher members than the singlet diffuse series. 
The 23P level is estimated to be overpopulated by a factor of 3. This could be caused 
by a dilution factor W = 0.05, which corresponds to R/r = 0.4, where R and r are 
the radii of the star and of the shell, respectively. The Si 1 triplet, 4p3P° — 4d3D, is 
found to be in emission in contrast with the well-known triplet 4s3S — 4p3P°, which is 
in absorption. Several strong lines are due to Fe ur. The effect of dilution in this 
spectrum is discussed. 


17 LEPORIS 


The spectrum of 17 Leporis presents a striking case of departures 
from thermodynamic equilibrium, presumably produced by dilution 
of the exciting stellar radiation in an extended atmosphere.’ The 
ordinary photographic spectrum has been described previously.? A 
spectrogram obtained with the Cassegrain spectrograph of the 
McDonald Observatory on March 5, 1939, covers the region A 3300— 
Ha. A reproduction of the violet part of this spectrogram is shown 
in Plate XVIIa. Ha is a strong bright line having absorption borders 
on both sides. This confirms the existence of a large nebulous shell 
around the star—which had already been inferred from a study of 
the absorption lines. 

Table 1 gives a list of those lines which are appreciably strength- 
ened in 17 Leporis, as compared with a Cygni. The wave lengths in 
the first column are approximate and have been taken from the list 
of lines of a Cygni by Wyse,’ whenever they were given in that 
paper. A few faint lines have been measured separately and com- 


* Contributions from the McDonald Observatory, University of Texas, No. 14. 


™Struve, Proc. Amer. Phil. Soc., 81, 241, 1939. 
2 Struve, Ap. J., 76, 85, 1932. 3 Lick Obs. Bull., No. 492, 1938. 
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TABLE 1 


LINES WHICH ARE ENHANCED IN 17 LEPORIS 


Most Impor- Int. 
Lower E.P 

tant Element Sun 
Ti 0.028 3 
Ti 0.012 5 
28038. 2... © .000 
S257 | Ti 0.028 5d? 
3301.5 Crit 2 
3394-3... Ti 0.012 3 
$403.4... Crit | 2.424 2 
3421.2 | Cri 2.411 4 
3442.0 Mnu 1.768 6 
3444.3... 0.150 4 
3450.6.. 2.052 
| 1.802 4d? 
3401.6.. Ti | 0.134 5 
3474.1... | Mnu | 1. 802 2 
3488.7 | Mnu 1.840 | 4 
3491.1 Tit 0.112 5 
3495.8 Mn it | 1.847 2 
3504.9 | 1.884 3 
3510.7 1.885 5 
1 V 3 
3520.1... Ti 2.039 2 
3550.9.. Vi 4 
S611 0.022 4-6 
| Fer 0.855 30 
Vi 1.067 5d? 
3599.0.. Ti 1 0.605 4 

| 2.695 IN 
5 

3613.7 SOE 0.022 4 
3618.8. . Fer 0.986 20 
3624.6.. Ti 1.216 5 
3630.8. . Se 0.008 4 
2. . 4 
3642.9f Scu ©.000 3 


Multiplet 


4, 


a4D,, —z4P%,, 
aiD,, 
aD,—z25P; 


b4F ,,, 
b?P,,,—y?P%,, 
asD,—25P3 
asD,—z5P3 


b4F,,, 
asD,—z5P3 
asD, —25P5 
asD,—25P5 


b?P,, —x2D3,, 


aF,—23D; 


ask; —25G8 


a‘P,, 

asF,—z5G3 


a?P,,—z?S2 


* Very strong. 
t Mostly due to Sc 11 in 17 Leporis. 
t Mostly Sc 11 in 17 Leporis but May be Cr in a Cygni. 
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TABLE 1—Continued 


tant Element Sun 

20567. Ser 0.008 4 a3D,—23F9 
3050.7... Tit 575 5 b?D.y,—y?F 3, 
3062.0. . 559 5 —y?F 3, 
3077.9.. Cru 2.693 a‘P,,, 
0.605 rod? a?F;,, 
Fet 0.051 asD,—25F3 

Call 6d? 578% 
3715.3 Cri 3.090 2 b4!D;,, —Z4F 
3720.2 Fet ©.000 40 asD,—25F3 

3737-9. Call 2.137 5 47Piy, — 
3741.5 Ti u 1.575 4 b?D.,,— y?D3, 
3745.7 Fet 0.087 8 asD,—25F3 
3759.4 Tiu 0.605 12d? 
3761.4. . Ti u 0.571 7 a?F,,, —2°F 
3900.6. 1.126 § a?Gyy, 


pared with new measurements of a Cygni made on McDonald Ob- 
servatory spectrograms. The identifications are also mostly those 
given by Wyse. The strongest contributing line has been listed in 
each case. In a few cases the relative importance of the blended 
components in 17 Leporis differs appreciably from that in a Cygni. 

The most significant fact revealed by Table 1 is the great strength- 
ening of 77 11. Apparently all multiplets are enhanced in 17 Leporis. 
All observable lines start from metastable levels, except a*F — z'G°, 
which starts from the ground level. The lines of this multiplet are 
probably enhanced even more than those of other multiplets which 
start from excited metastable levels. In comparing the lines of 17 
Leporis with those of a Cygni, it should be remembered that prac- 
tically all lines of the former fall on the 45° branch of the curve of 
growth.’ The curve of growth of a Cygni is not accurately known; 
but the turbulence is certainly less than in 17 Leporis, and only the 
weaker lines fall on the 45° branch.’ Because of the uncertainty in 
the curve of growth of a Cygni, it is not possible to ascertain definite- 


4Struve and Elvey, Ap. J., 79, 428, 1934. 5 Ibid., p. 427. 
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ly whether the low-level lines are enhanced more than those of higher 
excitation potential; but there is a suspicion that such may be the 
case. The first two ionization potentials of 77 are 6.8 volts and 
13.6 volts. 

We next notice the remarkable enhancement of the lines of Mn m1. 
All members of the multiplet a‘D — z5P° are greatly enhanced in 
17 Leporis. The lower level is metastable, and its excitation poten- 
tial is 1.8 volts. The first two ionization potentials of Mn are 7.4 
volts and 15.7 volts. 

Even more striking is the great enhancement of the line of Sc 11. 
All belong to the multiplets — and — both of 
which originate from the ground level. The line \ 4247 (a'D — z'D°), 
the lower level of which is metastable, with an excitation potential 
of o.3 volt, is also enhanced in 17 Leporis, although perhaps not 
quite as much as the lines of the multiplet aD — zF°. The first two 
ionization potentials of Sc are 6.6 volts and 12.8 volts. 

The lines of Cr 1 show a slight enhancement in 17 Leporis, but 
the effect is not very conspicuous. All lines originate from metastable 
levels. The first two ionization potentials are 6.7 volts and 16.6 
volts. There are some indications that V 1 lines are enhanced and 
that Nz 11 is about as strong as in a Cygni. These lines start from 
metastable levels. 

Of particular interest is Fer. Several lines of the multiplet 
a’D — 25F°, starting from the ground level, are greatly enhanced in 
17 Leporis. Appreciable enhancement is also observed in the multi- 
plet aSD — z'D°, originating from the ground level, and in a‘F — 
y°F°, originating from a low metastable level of excitation potential 
0.9 volt. The strongest lines originating from higher levels——for 
example, 4045 (a3F, — y*F%), excitation potential 1.5 volts—-are 
slightly weakened in 17 Leporis. But so are also the fainter mem- 
bers of the low-level multiplets—for example, » 3927.935 (a°D, — 
z°D$). Hence the effect is certainly in part produced by differences 
in the curve of growth. Whether there remains an excess in favor 
of the lowest multiplets is not certain. 

The lines of Mg 1 3838, 3832, and 3829 (3°P° — 33D) are of special 
interest in the study of departures from thermodynamic equilib- 
rium.’ The lower levels, 33P? and 3P%, are metastable, while level 


6 Struve, Proc. Amer. Phil. Soc., 81, 235, 1939. 
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35P? combines with the ground level of the singlet system, giving 
the intersystem line \ 4571.104 (3'S, — 3°P®). Hence, we should 
expect that the level 3°P? would be reduced in population when the 
radiation is diluted. The observational data for Nova Herculis were 
inconclusive,° but in 17 Leporis the three lines are clearly resolved 
from the hydrogen line H,. There is no appreciable weakening of 
d 3832. It is considerably weaker than \ 3838; but, on the other 
hand, it is stronger than \ 3829. Hence it is reasonable to attribute 
the slight differences from the relative intensities in a Cygni, where 
d 3832 is nearly equal to A 3838, to differences in the curves of 
growth. Evidently dilution does not sensibly affect the population 
of state 3°P%. 

If we consider three states—3'S,, 33P?, and a composite high level 
which includes the high discrete states—as well as the state of 
ionization, we have for the relative populations’ of states 33P? and 
3'So 

, +3 + ap; 


—=NH 12 
Ny y¥+3+ aWp,, 


where 
2 


R 
W = the dilution factor — , 
4r 


y= Ay 
Ax 
Pius kT 
—hv2; 
=e 


The lifetime of the state 3°P? in Mg1 is approximately 4 X 1075 


sec.* That of the combined upper state may be assumed to be of 
the order of 10~* sec. Hence we have, roughly, 


a=4X 


7 Proc. Amer. Phil. Soc., 81, 221, 1930. 
8 A.C. J. Mitchell and W. M. Zemansky, Resonance Radiation and Excited Atoms, 
p. 147, New York, 1934. 
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Assuming for the difference in excitation potential of levels two and 
three, 3 volts, and T = 15,000°, we have 


= 


In order that the population of level 33P,° be essentially that given 
by the Boltzmann formula, it is necessary that 


aW p23 >> + 3. 
The value of y is of the order of 1. Hence in our case 
W > 1074. 


This inequality is in good agreement with other estimates of WW. 
From the absence of strong emission lines produced by the expand- 
ing shell we infer that the radius of the shell cannot be excessively 
large. In other, similar stars an estimate of W = 0.01 was deduced 
from the helium lines.? 

It would be extremely interesting to investigate the intensities of 
the Mg 1 lines in shells for which W < 10~4, but at the present time 
no such observations are available. 

Table 2 lists those lines which are appreciably weaker in 17 Le- 
poris than in a Cygni. Most conspicuous is the multiplet of S71 
sp? 7D — 4’P°, AA 3854, 3856, and 3863. The ground level is odd. 
Hence, the lower level of the multiplet is normal, and the theory fully 
accounts for the weakness of the lines. Very faint traces of all three 
lines are present, and a value of W between o.1 and 0.01 seems fairly 
reasonable. The faintness of \ 4128 and d 4131, 37D — 4’F°, has 
been pointed out previously. The ionization potentials of Sz are 
8.12 volts and 16.27 volts. 

The only other lines in Table 2 are due to Fe 11. All but one of 
these arise from normal lower levels (c?F, d?D) which are higher than 
the lowest odd levels (z°D°, z*F°, etc.). These lines should all be 
faint. The great majority of Fe 1 lines in the ordinary photographic 
region arises from metastable levels. The strongest of these lines, 
d 4233, b4P,; — z4D%, is of about the same intensity in 17 Leporis 


9Struve and Wurm, 4p. J., 88, 105, 1938. 
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and in a Cygni; but the fainter lines are generally weaker in 17 
Leporis. The one anomalous line in Table 2, \ 3783.4, arises from a 
metastable level. The line is strong in a Cygni and, if the identifica- 
tion is correct, should be fairly strong in 17 Leporis. No satisfactory 
explanation has been found for this anomaly. The first two ioniza- 
tion potentials of Fe are 7.8 volts and 16.5 volts. 

The extreme weakness of the ultimate lines of Sv 11 4078 and 4216, 
5°S — 5?P°, has already been commented upon.’”® These lines are 
moderately strong in a Cygni, and in spite of differences in the curve 
of growth, we should have expected to find them reasonably strong 
in 17 Leporis. The ionization potentials are 5.6 volts and 1o.1 volts. 


TABLE 2 
LINES WEAKENED IN 17 LEPORIS 
d Lower E.P. Multiplet 
tant Element Sun 

ORO Fe 5.545 2 c?F;,,—x?Gyy, 
{Fett 5.930 abs. — x?F 3, 
/ 2% 3% 
| Fe 2. 766 ° bsP,, 


* This may not be all due to Fe 11. 


The ultimate lines of Ca 11 are about equally strong in 17 Leporis 
and a Cygni. The latter also shows the lines \ 3706 and A 3737, 
4’?P° — 57S, but in 17 Leporis both are heavily blended with Fe 1. 
This is unfortunate, since the Ca 1 lines would otherwise furnish an 
interesting test of the theory. 

Table 3 summarizes the results. The first part of the table con- 
tains those elements which are strengthened in 17 Leporis as com- 
pared with a Cygni, or which are of comparable intensity in both 
stars (Mgi, Fett metastable). The second part shows those ele- 
ments which are abnormally weak in 17 Leporis. The great weak- 
ness of Sr 11 cannot be explained by the theory of cycles. It is prob- 


1 Proc. Amer. Phil. Soc., 81, 242, 1939. 
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ably due to the fact that the second ionization potential is low, so 
that most atoms of strontium are doubly ionized. This, however, 
does not seem compatible with the great strengthening of 77 11, 
Sc u, and Fe 1 and with the weakening of Fe 11. However, the pres- 
ent theory is quite inadequate to explain the observational details. 
Perhaps the best procedure would be to determine the degree of 
ionization from the ultimate lines of Cat (\ 4226) and Cau (H 
and K). It will then be possible to study the relative abundances of 
those elements for which we observe the resonance lines (77 0, 
Sc u, Fe 1, and Sr 1) and then attempt to infer whether the intensi- 
ties of those lines of Feu which start from metastable levels are 
normal. But this work must be postponed until an accurate curve of 
growth for a Cygni is available. 
TABLE 3 
STRONG LINES IN 17 LEPORIS 


. 17 Leporis 
(in Order of Decreasing Ratio Joep) 


a Cygni 
Element Tin| Mnu Seu Cru Niu Vu Fer Mgt Feu 
724 6.6 6:7, 7-6 6.8 7.8 7.6 7.8 
Second I.P../13.6} 15.7 12.8 16.6 18.2 14.7 16.5 15.0 16.5 
Average E.P.| 1 1.8 | 0-0.3 2-3 I o-! 2.7 2-3 
WEAK LINES IN 17 LEPORIS 
: . 17 Leporis 
(In Order of Decreasing Ratio — = ) 
a Cygni , 
Feu 
Mg 
Element (Normal) ii Meu 
Average E.P....... 5-6 7-10 9 ° 


P CYGNI 
An excellent spectrogram of this star on Process emulsion was 
obtained with the ultraviolet Cassegrain spectrograph of the 
McDonald Observatory, on May 26, 1939. The linear dispersion 
varies from 21.7 A/mm at \ 3300 to 52.9 A/mm at A 4300. The star 
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TABLE 4 
ULTRAVIOLET LINES IN P CYGNI 
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(MEAS.) 


3451.62 
3455.00 
3459.32....--- 
3409.41..... 


to 


IDENTIFICATION 


Elem. 


Fe 111 


Fe 111 


Fe 


Fe 


Fe 111 


Fei 
Fe ut 


Het 


(Lab.) 


6 


6 


88 


08 


81 


22 


gI 


Transition 


;—25Gé 


as F, 


asF,—z5G; 
Unclassified 


Unclassified 


asF,—z5G;3 
Unclassified 


2$—7'P 
atF iy, 

Unclassified 
3p'So 


21S—6'P 


23P—213D 
23P — 203D 
23P—193D 
23P —183D 
23P—173D 


23P —163D 


R.V. 
(CORRECTED 
To Sun) 


* Interstellar. 


ABS. OR 

| 12) a 

| 4) 11 

3274.72...---- A 7 J —108 

3204 000 E 5) 

| E 2) 4 + 20 

| 2 | 9.80 + 70 

| 8 \ [—122 2 

A 4s Tiu 3.76 — 10 

3396.48.......| E 2 Feut | 6.71 = 

| 

A 4 Nu 7.16 — 89 

O08) A I Het 0.22 — 102 

A 2 Het 3.21 —122 _ i: 

i 

A 2 Hei 6.79 — 87 2 

| <A 2 Het 0.04 —124 

A Het — 184 

A 2 Het —I9gI 

| 

= 

= 

| 
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TABLE 4—Continued 
IDENTIFICATION 
d* (Meas.) Int 
Elem. dX (Lab.) Transition 

A 3 Het 8.97 23P —153D 
E 4 Sim 6.93 Unclassified 
3497.45....... A 5\ 3P— 133 
rf Het 8.64 23P —133D 
A Fu 3p5P1— 3d5Do 

E Fe 11 11.75 
SEOSUAOs A 3 Fi 5.61 3p5P;—3d5D3 
E 3 Fe 11 6.93 Unclassified 
A Het 0.49 23P—113D 
Het 4.50 23P —103D 
E 4 Het 23P—93D 
3589.79....... E 6 Si 0.46 4p'P,—4d'D, 
3603.76....... E 8) Fe 3.88 
E 6f | Her | 3-64 
Het | 4.28 23P —83D 


R.V. 
(CORRECTED 
TO Sun) 


— 94 


¢ He 3487.97 visible as absorption line on either side of the Si 111 emission line. 

t Probably an absorption line to violet of this emission line. 

§ He 1 3599.39 may be blended with this line. 


\ 
{ 
| 
| 


| 
— 
| 
\+ 
\ 
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| 
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: 
| 
f 
l+ 7 
Ios 
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: 
—119 
40 
130 
8 
f-117 
— 
4 


17 LEPORIS AND P CYGNI 737 
TABLE 4—Continued 
om IDENTIFICATION RV. 
(MEAS.) INT (CoRRECTED 
Elem d (Lab.) Transition To Sun) 
BOSC A 4 Het 2.06 
2662.68. .....- E 3 Fe 2.65 Unclassified |.......... 
3672.09. A 4 Ht 3.76 2-23 —12I 
3674.66....... A 4 Ht 6.36 2-22 —123 
9607 6 Ht 9.36 2-21 —124 
3081. 26. A 6 Ht 2.81 2-20 —I1I0 
A 8 Ht 6.83 2-19 — 131 
A 8 Ht 1:85 2-18 —125 
E 3) Ht 7-15 \- 8 
3703.14.. A 12 Het 5.07 2aP—7D ~=.......... 
3705.08. . E 9 Ht 3.86 
3722.09 E Ht (+ 28 
3748.46 A its 0.15 {—119 
3750.14....... E 35 
3768.86. . A II | i f—125 
3790.08. E 2 Si 11 1.41 4p3P.— 4d3D, — 18 


} 


“4 
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TABLE 4—Continued 


ABS. OR 


(MEAs.) Ew. INT 
BIO E 10f 
A 4 
3858.74....... A 5n 
3865 .35 A 4n 
3875.92. A I 
3881.10. A I 
3886.30. . A 25 
3924.92. A 6 
3033 A 8s 
A 2s 
3957-48. A I 
3902.58... A 12 
E I1f 
3967.92. A 25 | 
3970.02 E 25f 
3993 -59....... A 9 
3904.96....... 
A 7\ 
A 20| 
BOIS BS E 15f 


IDENTIFICATION 


Elem. 


Ht 


St ul 


Het 


d (Lab.) 


7.90 
6.56 


9.61 


8 


Transition 


2-10 
4p3P.— 4d3D, 


23P—63D 


23P—63S 

2'P—g'D 

e4D;,, —z4D;,, 


23S — 33P 
2-8 


2'P—8'D 
47Sy, % 
383P,—3p'P: 


2-7 
3s'P,—3p'P; 


2'7P?—7'D, 


23P — 53D 


R.V. 
(CORRECTED 
TO Sun) 


f—116 


|| Interstellar. 


| 
I 5 39 2 9 
ro 
Het 7.54 154 | 
H 8 | 
Het 8.18 —159 
Ou 2.19 — 76 
/ 
Ht 9.03 | 
Het 6.53 107 | 
3-97 4 | 
Vu 5.85 —82 
| 
4-73 + 6 
{—147 
f— go 
Het 2 | 
—132 
| 
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spectrum is measurable from \ 3265 to Hf. Since the region from 
HB to about 3900 was studied in detail by Struve™ on a plate of 
somewhat greater dispersion, the present measurements were limited 
TABLE 5 
He LINES IN P CYGNI 


| 
| ABSORPTION EMISSION ABSORPTION EMISSION 
\ n n 
Int. | Vel. Int. | Vel. | Int.| Vel. Int. | Vel. 
2'S—n'P 23P 
5015 3 7 | 160 | — 6 4713 4 7 | —136| 10 | — 7 
3004 4 |12 | —147] 11 | + 6 4120 5 & | 5 — 2 
3613 5 8 | —123] 6] +17 3867 6 40) 
3447 6] 4 | —125 7 | Masked by H 
3354 7 4nj —122) I — 2 3734 
3652 8 | 3 
2'P—n'S 3509 9 | Masked by Fe III 
5047 4 I = 570 3000.95 
4168 6 1 | — 83 23P —n3D 
39035 8 447! | —158] 20 | — 8 
38109 6| 1 | —114) | + 5 
2'P—n'D 3705 7 | Blend with H 3704 
4921 4 | 10 | —154| 10 | + 6 3634 8 8 — o98| 6 | +16 
4387 5 | 10 | —123]} 6 | —12 3857 9 | Masked 4/+ 8 
4143. 6 g | —115| 2] — 8 3554 10 5 | —105| 1] +10 
4009 7 8 | —120] 2] — 6 3530 II 6 | —106} 2] +7 
| 3926 8 6 | —107] —I5 12 6 | — 92] 1| + 6 
3871. adh 3498 13 5 | — 86] 1 | +71 
3833 10 | Blend with H 3487 14 | Masked by Sz 111 
3805 11 | Just visible. May be 3478 15 
blend with S7 III 3471 16 2 | —I191 
| 3784 12 | This line is not 3465 17 2 | —184 
| 3768 13 | Blend with 7 3456 19 2| — 87 
3756 14 | Absent 3453 20 | 2] —122 
3450 21 I | —102 
23S —n3P 
3888 3; Blend with 
3889 


on the long wave-length side by He 4026. The results are assembled 


in Table 4 


Helium.—The lines of six He 1 series represented in P Cygni” are 
assembled in Table 5. A special search was made for the higher 


1 A p. J., 81, 66, 1935. 
'2 The data for the lines of longer wave lengths than \ 4026 are taken from Struve, ibid. 
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members of the two series, 2"P — n'D and 23P — nD. The exten- 
sion of the triplet series to m = 21 is a striking feature of the spec- 
trum. The singlet series extends tom S 11. These facts suggest that 
the 23P and the 2'P levels have relative populations which are not 
in accordance with their Boltzmann factors or, stated in another 
way, that there is a departure from thermodynamic equilibrium. 
The two levels are separated by 2048 cm~', which means that the 
relative Boltzmann factor, 2'P:25P, should be 1.2:1. For practical 
purposes we may assume the two levels should be populated in the 
ratios of their statistical weights, 3:1, if thermodynamic equilibrium 
holds. To test this point a comparison was made with the stars 
55 Cygni and + Scorpii. The results are shown in Table 6 and in 
Figure 1. 
TABLE 6 


LAST VISIBLE MEMBER OF THE 
DIFFUSE He I SERIES 


Singlet Triplet 


In the latter we have plotted against serial number n, the 
quantity log X%, which has been used by Goldberg'} and which is 
proportional to the oscillator strength, f. The broken lines for 
t Scorpii and 55 Cygni connect the last members of the 2'P — n'D 
and the 25P — n3D series. These lines are roughly parallel. Hence, 
if we start from the fact that in P Cygni the last singlet line is 
n = 11, we should, by analogy, expect that the last triplet would be 
n = 13.5. In reality we observe n = 21. The difference in the ordi- 
nates of the triplets m = 13.5 and m = 21 gives us A log NX = 0.58, 
from which we infer that the ratio of the populations 23P/2'P is 
four times larger than the value found for thermodynamic equi- 
librium. 

This result is complicated to some extent by the fact that Gold- 
berg has included in his expression for log X( the factor (A’/A)3, 
which allows for a yradual change in the general opacity of the 


I 


13 Ap. J., 89, 630, 1939. 
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atmosphere with . It is doubtful whether this law should be used 
for P Cygni. It is better to assume that the expanding shell is 
essentially transparent outside the lines. Since we are comparing 


log X 


n 
Fic. 1.—The upper curve represents the Het triplets; the lower curve, the He1 
singlets. 
very weak lines, we shall suppose that the conditions for absorption 
in a thin slab are fulfilled. The equivalent widths of the last lines 
of the two series are 


= Nox N tH f(2—21) 


mc? 


From the observations it follows that 


Hence, 
= 
or 
N, 


4 
q 

| 

+2 
| 
A 
+1 
° 4 
1, 4co 
4 oe 
“Ax 
55 “CYC 
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The values of f have been computed by Goldberg" only to 2 = to. 


Hence, we adopt the asymptotic expressions: 


f (triplets) = 4.69n"3, 


f (singlets) = 3.33273. 


The result is V,/N; = 0.17, while the thermodynamic value is 0.3. 
The overpopulation of the 23P level corresponds to a factor of 2. 

Under the circumstances the best that we can do is to take the 
mean resulting from the two methods used, giving an overpopula- 
tion by a factor of 3. This result shows fairly conclusively that the 
departure from thermodynamic equilibrium is present but that it 
is not excessively large. The computations by Struve and Wurm‘ 
show that for a temperature of 25,000° K the relative populations 
of the 23P and the 2'P states vary with the dilution factor W as 
follows: 


Hence the best value for the dilution factor of the He shell in P 
Cygni is 


W =0.05. 


This corresponds to a value of 


where R and ¢ are the radii of the star and shell, respectively. This 
unexpectedly small extent of the expanding shell is in good agree- 
ment with other evidence."® 


4 Ap. J., 90, 414, 1930. 
5 Ap. J., 88, 98, 1938. 


1 Struve, Proc. Amer. Phil. Soc., 81, 242, 1939. 
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The stronger absorption lines of Het cannot be used for the 
evaluation of V, and NV, because we have no knowledge of the curve 
of growth. In Figure 2 we have plotted the observed values of 
log A/d against Goldberg’s values of log X%. The systematic differ- 
ence between triplets (open circles) and singlets (closed circles) is not 
consistent with a single curve of growth for both series. The presence 
of strong emission lines in the triplets may account, in part, for this 
phenomenon. An inspection of the plate suggests that the singlet 


“3 


log A/X 


| 
—I +1 


log 


WwW 


Fic. 2——The diffuse triplets of 7e 1 are shown as open circles, the diffuse singlets 


as closed circles. 


absorption lines are systematically broader than the triplets. The 
diagram is directly comparable to those published by Goldberg for 
various stars. He did not observe the ultraviolet triplets, and it 
should be of interest to investigate whether normal B stars show 
the same anomaly. 

We have measured the equivalent widths of the He 1 absorption 
lines and the total emissions, referred to the energy included in one 
angstrom unit of the continuous spectrum, in the immediate vicinity 
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of the line. These values are given in the third and fourth columns | 
of Table 7. 
TABLE 7 


He 1 LINES IN P CYGNI 


Emissi (1/W)N | 
Transition E.W. (X1073) ge—x/k 
(X107}8) | 
23P— 43D...... 4471 0.815 2.3 1.4 
4026 .750 0.645 20 2.3 
| 3819 ©. 245 1.8 
3634 . 240 0.13 
108D........ 3554 | 
123D... 3512 101 0.063 5.8 5.4 | 
3408 
4388 730 0.012 0.07 Ours 
4144 440 0.011 0.13 
71)... 4009 
SD)... 3926 305 
3870 340 | 
23P— 538..... 4120 120 0.26 | | 
3013 0.135 0.09 | | 
Suppose that the expanding shell has a thickness /J and an aver- | 


age radius r. Then the total energy, in ergs, radiated by the entire | 
shell in one second, in all directions, is | 


We have measured this radiation in terms of the radiation of the 
continuous star spectrum. This is 


tRcp,dv , 


where | 
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We choose for dv the frequency interval 6y, corresponding to 1 A 
at the wave length considered. The quantity listed in the fourth 
column of Table 7 may be designated as x. Then 


nH hve Ana = bv . 


Recalling that, approximately, 


R? 
W=—, 
4r 
we have 
I p,6v 
— HN, = —— 
Ae 


The values of A,. are given by Goldberg.'* Hence we may compute 
(1/W)HN,,; these quantities are given in the fifth column of Table 7, 
where we have used T = 20,000°. This choice must be explained. 
The energy from the nebula and from the unit in which we have 
measured it is weakened in the same proportion by interstellar ab- 
sorption. Hence, in converting our unit into ergs per second, in all 
directions, we must use the amount of radiation within 6y as seen 
from a representative point within the nebula. Sherman and Mor- 
gan’ have shown that ‘“‘a considerable part of the observed redden- 
ing is probably due to physical conditions in the immediate vicinity 
of the star.” It is impossible to ascertain whether the light of the 
continuous spectrum is already weakened when it reaches the 
nebula. 

If it is, we should use T < 20,000°. This would greatly affect the 
unit but would not materially change the Boltzmann corrections. 
We should remember that for 7 = 10,000° our unit would be almost 
ten times smaller and that all values of the last two columns of 
Tables 7 and 8 would be ten times larger. 

The weaker absorption lines give the following quantity: 


- 


A, = 
mc? 


Ap. J., 89, 516, 1939. 
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We may now compare the values of NV, derived from the emission 
lines with those of V, derived from the absorption lines. We may 
write 


Nn _ 


N, g.b, 


where b, is Menzel’s'® factor, which measures the departure from 
thermodynamic equilibrium. This quantity is not known for helium. 
The last column of Table 7 gives the values of 


I 
W HN, 


gre 


These values should be compared with the corresponding quan- 
tity obtained from the weak absorption lines: 


Consider, for the diffuse triplets, the line \ 3554. The equivalent 
width is A, = 0.2 A. Hence, 


= 4 X 105, 


and 
N.H 


£2€ kT 


= to’. 


Table 7 gives the results obtained from the emission lines. 
Since the capture cross-sections for the various series of He I are 
not known, and since we do not even know to what extent recom- 


8 Ap. J .;88, 60, 1938. 
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747 


bination is responsible for filling up the various levels, we cannot de- 
termine 0,/b,. The tabular values show that the singlet D levels are 
relatively less populated than the triplet D levels; and since the 
differences in the statistical weights have already been taken into 
account, we must attribute the effect to differences in the values 


of b,. 


Hydrogen.—The Balmer series is followed in absorption to 
d 3073.76 (2-23) and in emission to \ 3697.15 (2-17). Microphotom- 
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Fic. 3.—Radial velocities of the absorption lines of hydrogen 


@ Emission lines | © Absorption lines 


Intensity of emission lines 


eter tracings strongly suggest that the emission lines also persist 
to the twenty-third member, but the eye does not distinguish the 
emission wings after number 17. In Figure 3 the apparent radial 
velocity of the absorption lines of hydrogen is plotted against the 
serial number of the line. The parallel plot of the emission line in- 
tensities indicates that the change in apparent velocity is caused by 
the relatively greater encroachment of the emission line on the 
violet absorption line for the early members of the series. This 
causes the absorption lines to be measured with systematically 
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shorter wave lengths and thus leads to a higher apparent negative 
velocity. 
The persistence of the Balmer series to number 23 is evidence of 
high luminosity, as pointed out by Unsold.’? From Unsdéld’s calibra- 
tion we find for the absolute magnitude about —5. Ina recent paper 
Sherman and Morgan"? found for P Cygni Myi, = —4.4. | 


TABLE 8 
BALMER SERIES INTENSITIES 


(1/W)N,,H 
= EW. Emission (1/W)N nll 

Intensity (Xt0-'5) 

E (10-16) 
5. °.820 3.88 0.8 3.0 
6. 0.915 1.33 0.8 2.2 
1.04 1.4 

8. 1.100 2:92 

Q. 0.775 2.4 
IO. ©.510 0.64 6.0 

II 0.440 0.20 3.0 | 
0.400 0.10 2.4 
0.390 0.13 4.0 4 
14. 0.200 0.16 8.4 
15 ©. 200 0.16 2.0 6.6 
16. 0.360 0.28 29.4 14.4 
0.140 0.17 24.0 10.4 
18. 0.060 ro 320 9.3 

1Q. 0.070 0.18 44.5 5 | 

20... 0.020 0.06 | 
0.030 0.10 41.1 2 
22. 0.010 0.03 15.6 4.1 
0.030 0.01 6.5 1.6 


* Blended with //e 1 3888.64. 


The measurement of the intensity of an absorption line flanked 
by an emission line, or vice versa, is complicated by the distortion 
of one by the other. The values of the equivalent widths of the ab- 
sorption lines and. the intensities of the emission lines, recorded in 
Table 8, second and third columns, were obtained by the measure- | 
ment of the undistorted half of each line. We assumed that the true | 
line is symmetrical. The fourth column gives, for the emission lines, 
the quantity (1/W)N,H, and the fifth column gives the quantity 


"9 Unpublished. 
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[((1/W)N,H/2n?e-”’*"]. There is a tendency for this quantity to in- 
crease from n = 5 ton = 19. The decrease from u = 19 ton = 23 
may not be real, because of the uncertainty in drawing the continu- 
ous spectrum and because of the weakness of the lines. The increase 
from 2 = 5 tom = 19 may be considered as an effect of the factor b, 
of Menzel. From the observations 0,,/b; = 5. The tables by Baker 
and Menzel'® show that for 7, = 20,000° neither their hypothesis 
A nor their hypothesis B provides for a ratio greater than about 3. A 
lower electron temperature, of the order of 7, = 10,000°, would be 
more consistent with the observations. It is premature to discuss 
this matter further. 


TABLE 9 
N.H FOR HYDROGEN 
n (Upper) E.W. log NHf log f* log NH NH 
20; 0.020 — 3.35 14.57 3.7X10%4 
.030 I1.40 14.72 5.2 
.O10 10.90 3.48 14.40 


* From Menzel and Pekeris, M.N., 96, 82, 1935. 


Table 9 gives the results for V,H from four weak H absorption 
lines. The mean is N,H = 5 X 10°; N,H/ane~”/*? = 2.3 X 10°. 
This is almost identical with the value of (1/W)N,H/2n?e~”** 
found for the emission lines having small values of m. Thus, the 
mean from Table 8 of the emission lines nm = 5,” = 6, and ” = 7 is 
2.8 X 10. Hence, 


1 bs 2.8 X 10% 
— — = ————— =] .2., 


W b, 2.3 X 10° 


From the tables of Menzel and Baker" we infer that ./b, is sensitive 
to T., as well as to the mechanism involved; but probably a value 
of be/b, = 3 is not excessive. This would give 


and this is, of course, not consistent with the limitation W < 1. 
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We conclude that the emission lines are in reality too weak to be 
consistent with the absorption lines. Several explanations are pos- 
sible: 

a) The unit of measurement of the emission lines may be in error 
by as much as a factor of 10. If we adopt T = 10,000° instead of 
T = 20,000°, we obtain, roughly, 


so that, for b6/b. = 3, we find W = 0.3. 

b) The assumption of a uniform shell, without appreciable ob- 
scuration by the star and without overlapping of emission and ab- 
sorption, may not be permissible. This assumption rested upon the 
evidence from Hei, C u, and Mg u that W = 0.05. 

c) The metastability of the 2s level of H may cause this level to 
be overpopulated. We are unable to test this hypothesis because 
we do not know enough about the lifetime of the 2s level. 

It is clear that the difficulty could best be overcome by observing 
the Paschen absorption lines in P Cygni. Attention must also be 
called to the fact that, although the inconsistency between absorp- 
tion and emission may be overcome in P Cygni, it will be even more 
serious in 17 Leporis. The emission lines in that star are very faint. 
They are visible only in Ha and Hf. The absorption lines, on the 
other hand, are almost as strong as in P Cygni. The dilution factors 
are fairly well ascertained in both stars. 

N II.—The nitrogen lines present an interesting case of de- 
partures from thermodynamic equilibrium: \ 3437 appears in ab- 
sorption only, while the entire multiplet AX 4601-4643 shows strong 
emission and absorption. 

N IIT.— 3745.83 may be present, but the absence of \ 3754 and 
d 3771 makes the identification doubtful. 

O II.—Present: \ 3727 (blended); \ 3882 (1); \ 3470 (blended). 

O IIT.—> 3961 is absent. 

F II.—The multiplet 3p°P — 3d°D° may be present, but two 
lines are blended with Fe ur. The radial velocity of the third line 


is abnormal. 
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Si III.—One of the most striking features of the spectrum of 
P Cygni is the peculiar behavior of the Si 1 lines. The triplet 
dA 4552, 4567, and 4574 is in absorption with an emission wing on 
the red side of \ 4552, while the lines in the ultraviolet are pure emis- 
sion lines with no evidence of any absorption. Table 10 summarizes 
the results. 

TABLE 10 


Si III IN P CYGNI 


d (Lab.) EorA Intensity Velocity Transition 
AT. E 2 18 4p3P.—4d3D, 
3700.12... Blended with Hip 4p3P,—4d3D, 
3806.56*. .. E 6 19 4p3P,.—4d3D, 
A+E 8,1 95 (—7) 4838S; — 4p3P, 
BAO A 7 103 483S;— 4p3P, 


* Attention was called to the peculiar intensity of this line by Struve (A p. J., 81, 95, 1935). 


The Grotrian diagram for Si 1m (Fig. 4) indicates the relationship 
between the absorption triplet and the emission triplet. The 43P 
level is common to both. The interpretation of this departure from 
equilibrium may be possible when spectrograms are available for 
shorter wave lengths, to include the multiplets 3d*D — 4p3P° 
(A 3086, etc.) and 4p3P° — 5s3S (A 3230, etc.) as well as the singlet 
transition 4p'P® — 5s'S (A 3185). We note the fact that the 3P° 
level can reach the lowest triplet level only by two steps, a fact 
which may have something to do with the anomalous behavior of 
the two triplets. 

S ITT.—> 3324.85 is present. 3338.32 is blended with N 11. 

Fe ITI.—The Fe 111 lines are a marked feature of the spectrum. An 
examination of Plates XVIIb and XVIIIa shows quite clearly that 
the great strength of the Fe 11 lines is not a characteristic of the 
early B spectral classes, as these lines are not visible in 7 Scorpii 
and 55 Cygni. In P Cygni they are, on the contrary, the strongest 
lines in the spectrum, with the exception of those due to hydrogen 
and helium. 

The presence of Fe 111 in P Cygni was pointed out by Swings and 
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Struve.”° The present study strengthens the identification. The 
Fe 11 spectrum has been analyzed by Swings and Edlén.”' 

In his paper on the photographic region of the spectrum of P 
Cygni, Struve™ measured a strong line at \ 4396, which he tentative- 
ly identified as the O m1 line \ 4395.95. It is now much more probable 


's 'p 'D 45 3p 3D bd 2 36 | 
o t 2 a 23 4 


100,000 


Excitation potentials in volts 


250,000 


Fic. 4.—Si mt: Energy-levels and spectral lines of astrophysical interest (from 
Bowen, Phys. Rev., 39, 8, 1932). 


that the star line is identical with Fe 111 4395.78. This removes the 
necessity of invoking abnormal fluorescence to explain the abnormal 
intensity of the O 11 line. 

The mean radial velocity from the absorption lines of Fe 1 is 
—123 km/sec. This value agrees remarkably well with the correla- 


20 Ap. J., 88, 621, 1938. 
21 Dr. Swings has very kindly placed at our disposal, in advance of publication, the 
results of the analysis of Ke 1 by himself and Edlén. 
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tion found by Struve'! between ionization potential and radial 
velocity. With its ionization potential of 30 volts Fe 11 should be 
similar to Het and Ni. The ratio of emission to absorption for 
the Fe m1 lines agrees with the relation between this quantity and 
the ionization potential.” 

Swings has called attention to the remarkable difference in the 
intensities of Fett 4420 and 4164. These lines are almost 
identical in y Pegasi, while in P Cygni \ 4420 is very much stronger. 
The latter starts from a metastable level, while \ 4164 starts from a 
normal level. We estimate that W < o.1. 

Although the strong lines of Fe ur in the ultraviolet regions origi- 
nate from levels which are not strictly metastable, it is clear that 
they must behave very much as though they were metastable. These 
lines originate from the levels aSF and aF. The former can combine 
only with a single lower odd level, 4p’P°. The line would be very 
weak, because of the change in multiplicity and the large change in /. 
Level a’F can combine with two lower odd levels, 4p7P° and 4p°P°. 
These intersystem combinations would also be weak. Hence the 
Fe 111 lines must be very strong for the same reason that \ 3832 of 
Mg tis strong in P Cygni. 

Interstellar lines.-Two strong absorption lines are of interstellar 
origin: 3383.76 (77m) and 3933.67 (Ca 1). 

TABLE 11 


UNIDENTIFIED LINES IN P CYGNI 


» (Star) EorA Intensity d (Star) EorA Intensity 
3405.50 E 2 ROC A Is 
3574.08 A I |: A I 
2047 ... A 6 A 3 
3049.29 E 6 


YERKES OBSERVATORY 
AND 
McDonaLp OBSERVATORY 
August 1939 


2 Ab. J., 81; 87, 1935. 
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Molekiilspektren und Molekiilstruktur, Part 1: Zweiatomige Molekiile.. By 
Dr. GERHARD HERZBERG. Dresden and Leipzig: Theodor Steinkopff, 
1939. Pp. xv+4o04. Figs. 169. Tables 37. Rm. 21; bound, Rm. 22.50. 


This book concerns itself only with diatomic molecules. Part II will 
be devoted to polyatomic molecules and will appear as soon as possible. 
Part I is consistent in plan with Dr. Herzberg’s earlier and well-received 
book, Atomic Spectra and Atomic Structure, and it should meet with equal 
favor. 

The author has attained his objective of providing a suitable introduc- 
tion for the beginner as well as a helpful compilation for the research work- 
er. The style is lucid, and the observational material has been brought up 
to date. Both the empirical and the theoretical aspects of the problem are 
treated. In those instances where proofs are deemed too difficult for the 
beginner, an attempt is nevertheless made to provide him with an insight 
into the situation. He receives further assistance from the use of fine 
print in which are given the details of those portions of the book which 
may be omitted in an initial survey of the field. Good diagrams, figures, 
and tables are abundantly employed. 

After a review of the fundamentals of atomic theory, an excellent 
presentation is made of the empirical relationships inherent in diatomic 
molecular spectra, covering the entire spectrum from the ultraviolet 
through the infrared. This is followed first by an adequate analysis of the 
rotation-vibration problem, made pointed by the use of examples from the 
literature. Electronic states and transitions are examined in the two suc- 
ceeding chapters. The remaining chapters are devoted mainly to a dis- 
cussion of the electronic configurations possible on theoretical grounds and 


« [While this review went through the press, there appeared a revised translation of 
it under the title Molecular Spectra and Molecular Structure (Prentice-Hall, Inc., 1939; 
price, $6.50). The translation was made with the co-operation of the author by Pro- 
fessor J. W. T. Spinks of the University of Saskatchewan. Several additions and 
changes were made for the translation, and the table of molecular constants for the 
ground states of all diatomic molecules thus far investigated has been brought up to 
date as of July, 1939.—EbITorS. | 
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the bearing of this analysis upon heats of dissociation, dissociation proc- 
esses and spectra, and valence. The book concludes with a chapter of ex- 
amples, observational results, and applications to the fields of physics, 
chemistry, and astrophysics. 

| References to the literature are thorough, and the indexes are usefully 


done. 
ARTHUR ADEL 


Lowell Observatory 
Flagstaff, Arizona 
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